SOUND  PHASE  AND  AMPLITUDE  FLUCTUATIONS 
IN  AN  ANISOTROPIC  OCEAN 


Charles  Homer  Alexander 


Library 

Naval 

Monterey,  Caftfcmii  939WT 


NAVAL  POSTGRADUATE  SCHOOL 

flonterey,  California 


T 


SOUND  PHASE  AND  AMPLITUDE  FLUCTUATIONS 
IN  AN  ANISOTROPIC  OCEAN 

by 
Charles  Homer  Alexander 


Thesis  Advisor: 


H.    Medwin 


December    1972 


T 


Appiovnd  foK  pubtic  fintzaAz;   dihtnlbuXiovi  tmtani£e.d . 


i 


Sound  Phase  and  Amplitude  Fluctuations 
in  an  Anisotropic  Ocean 


by 


Charles  Homer  Alexander 
Lieutenant,  United  States  Navy 
B.S.,  Purdue  University,  1964 


Submitted  in  partial  fulfillment  of  the 
requirements  of  the  degree  of 


MASTER  OF  SCIENCE  ENGINEERING  ACOUSTICS 

from  the 
NAVAL  POSTGRADUATE  SCHOOL 
December  1972 


A  37*3 
c-' 


Library 

Naval  Postgraduate  School 

Monterey,  California  93940 


ABSTRACT 

Sound  of  constant  amplitude  and  frequency  was  trans- 
mitted simultaneously  in  three  orthogonal  beams  over  a 
distance  of  1.5  meters  in  the  upper  ocean.   Time  records 
of  the  resulting  phase  and  amplitude  fluctuations  of  the 
sound  beams  were  studied  by  means  of  auto  and  cross  spectral 
analysis  and  correlation.   The  time  lag  between  correspond- 
ing peaks  of  the  phase  fluctuation  autocorrelation  functions 
of  vertical  and  horizontal  beams  indicate  movement  of  in- 
homogeneities  between  the  sound  fields  due  to  water  particle 
motion  caused  by  surface  wave  action.   Envelope  correlation 
times  of  the  phase  fluctuations  are  found  to  be  approximately 
one-half  as  great  in  the  mixed  layer  as  in  the  thermocline, 
and  are  greater  in  the  vertical  than  in  the  horizontal 
direction  in  the  thermocline.   Anisotropy  in  the  thermo- 
cline is  also  indicated  by  the  variance  of  phase  fluctuation 
being  greater  for  sound  paths  in  the  horizontal  than  in  the 
vertical  direction.   The  autocorrelation  functions  of 
amplitude  and  phase  fluctuations  in  any  one  direction  are 
similar. 
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I.   INTRODUCTION 

A.   BACKGROUND  AND  OBJECTIVES 

When  a  sound  wave  is  propogated  in  the  upper  layers  of 
the  ocean,  both  its  amplitude  and  phase  fluctuate.   These 
fluctuations  result  from  the  motion  of  the  medium  in- 
homogeneities  in  the  sound  path,  and  motion  of  the  medium 
itself.   Relatively  little  study  has  been  done  in  the 
relationship  between  the  ocean  microstructure  and  sonic 
fluctuations.   This  situation  exists  because  early  re- 
search on  the  interaction  of  sound  and  the  ocean  medium 
concentrated  on  the  effect  of  the  ocean  medium  on  the  sound 
wave's  intensity  only;  it  wasn't  felt  necessary  to  inquire 
about  how  the  sound  beam's  amplitude  and  phase  responded 
to  the  medium.   However,  a  knowledge  of  how  the  ocean  medium 
effects  the  phase  and  amplitude  of  a  sound  beam  would  now 
be  useful  from  two  standpoints.   Devices  are  being  developed 
for  use  in  the  ocean  that  utilize  sound  amplitude  and  phase 
information  (acoustic  imaging  devices,  underwater  communica- 
tion devices  etc.)   A  knowledge  of  the  effects  of  the  ocean 
medium  on  the  sound  amplitude  and  phase  is  therefore 
mandatory.   A  knowledge  of  the  relationship  between  the 
nature  of  the  ocean's  microstructure  and  the  amplitude  and 
phase  fluctuation  of  a  sound  beam  passing  through  that  med- 
ium could  enable  one  to  derive  information  about  the  ocean 
microstructure  by  measuring  such  fluctuations. 
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The  objective  of  this  thesis  was  to  take  in-situ 
acoustical  measurements  over  a  range  of  approximately  1.5 
meters  in  the  near  surface  region  of  the  ocean  in  order  to 
study  the  relationship  between  sound  phase  and  amplitude 
fluctuations  and  the  nature  of  the  medium.   To  do  this,  a 
constant  frequency,  constant  amplitude  sound  beam  was 
transmitted  in  each  of  three  orthogonal  directions 
simultaneously.   Measurements  of  the  position  of  the  thermo- 
cline,  swell  magnitude  and  direction,  and  other  environmental 
parameters  were  made.   Spectral  analysis  and  cross-spectral 
analysis  are  performed  on  the  resulting  phase  and  amplitude 
fluctuation  records,  and  the  sound  phase  and  amplitude 
fluctuations  are  related  to  the  behavior  of  the  medium. 

B.   CAUSES  OF  FLUCTUATIONS 

When  a  continuous  sound  beam  of  constant  frequency  and 
constant  amplitude  is  transmitted  through  a  volume  of  ocean, 
it  is  found  that  the  sound  beam  arriving  at  a  receiver 
placed  some  distance  from  the  transmitter  will  have  a 
fluctuating  amplitude  and  a  fluctuating  phase.   The  causes 
of  these  fluctuations  are  described  below. 

1.   Sound  Amplitude  Fluctuations 

A  sound  beam  traveling  through  an  ocean  medium 
encounters  temperature  and  salinity  inhomogeneities  and 
"patches"  of  small  bubbles  of  various  sizes.   The  in- 
homogeneities and  patches  vary  spatially  with  the  motion 
of  the  medium    (Medwin,  1970)  .  The  inhomogeneities  present 
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a  change  in  acoustic  impedance,  rho»C,  to  the  sound  beam, 
and  the  beam  is  partially  reflected  from  the  inhomogeneity, 
and  partially  transmitted  through  the  inhomogeneity  at 
some  refraction  angle.   If  a  sound  beam  traveled  through 
an  ocean  volume  in  which  the  inhomogeneities  and  bubble 
patches  were  "frozen"  in  place,  reflection  and  refraction 
would  still  take  place,  but  the  received  sound  beam  would 
have  a  constant  amplitude.   However,  when  the  "frozen" 
inhomogeneities  are  allowed  to  move  about,  to  move  into 
and  out  of  the  sound  beam's  path,  then  the  reflection  and 
refraction  of  the  sound  beam  changes  with  time,  and  hence 
the  received  sound  amplitude  will  fluctuate  accordingly. 
The  motion  of  temperature  and  salinity  inhomogeneities  and 
bubble  patches  into  and  out  of  the  sound  beam's  path  causes 
sound  amplitude  fluctuation,  due  to  the  sound  beam  being 
both  reflected  and  refracted  by  the  changing  acoustic  im- 
pedance of  the  medium.   The  extent  of  the  fluctuations 
depend  upon  the  impedance  mismatch  between  the  medium  and 
the  inhomogeneities,  size  and  shape  of  the  inhomogeneities 
or  patches,  and  the  frequency  with  which  they  move  into 
and  out  of  the  beam's  path. 

The  bubbles  mentioned  above  can  have  another  effect 
on  the  sound  beam.   If  the  moving  bubbles  have  a  resonant 
frequency  at  or  near  the  frequency  of  the  sound  beam,  then 
the  bubbles  will  resonate,  absorb  and  scatter  the  incident 
sound  beam.   This  reduces  the  sound  beam's  intensity,  and 
causes  the  amplitude  to  fluctuate  with  time. 
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2 .  Sound  Phase  Fluctuation 

A  fluctuation  in  phase  of  the  received  sound  beam 
results  from  a  change  in  the  integrated  speed  of  sound 
propagation,  C,  between  the  source  and  the  receiver.   Such 
a  change  in  C  in  the  ocean  results  from  the  following: 

a.  A  change  in  temperature,  salinity,  and/or  the 
resonant  bubble  population  in  the  volume  of  water  through 
which  the  sound  beam  is  traveling  will  effect  the  speed  C. 
This  will  occur  when  temperature  or  salinity  inhomogeneities 
and/or  resonant  bubbles  move  into  or  out  of  the  path  of 

the  sound  beam. 

b.  Motion  of  the  medium  with  a  vector  component 
in  the  direction  of  sound  propagation  will  add  to  or  sub- 
tract from  C    (doppler  shift) .   The  entire  ocean  is  in 
constant  motion;  the  predominant  motion  encountered  in 

the  surface  layers  of  the  ocean  is  generally  orbital  motion 
due  to  surface  wave  action. 

3.  Sound  Amplitude  and  Phase  Fluctuation 

The  cause  and  effect  relationship  of  the  ocean 
medium  and  sound  amplitude  and  phase  fluctuations  can  be 
summarized  as  follows: 

If  a  sound  source  and  receiver  are  positioned  in 
the  upper  layer  of  the  ocean,  the  orbital  motion  of  the 
volume  of  water  between  the  source  and  receiver  (in  response 
to  the  surface  wave  action)  carries  temperature  and  salin- 
ity inhomogeneities  and  varying  bubble  populations  into 
and  our  of  the  sound  beam's  path.   Amplitude  fluctuation 
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results  due  to  scattering  and  refraction  of  the  sound  beam 
as  it  encounters  the  differing  acoustic  impedance  of  the 
inhomogeneities  and  bubble  patches .   Phase  fluctuation 
occurs  due  to  a  change  in  the  integrated  speed  of  sound 
between  the  source  and  the  receiver,  and  also  due  to  the 
medium  velocity  component  in  the  direction  of  sound 
propagation. 
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II.   THE  EXPERIMENT 

A.  OBJECTIVE  OF  THE  EXPERIMENT 

The  objective  of  the  experiment  was  to  propagate 
constant  amplitude,  constant  frequency  acoustic  signals 
in  three  orthogonal  directions  simultaneously  at  various 
depths  in  the  upper  layer  of  the  ocean,  and  to  record  the 
resulting  amplitude  and  phase  fluctuations.   Frequencies 
used  in  the  experiment  were  65  kHZ,  and  105  kHZ. 

B.  INSTRUMENTATION 
1.   The  Frame 

a.  Construction 

A  rigid  frame  was  constructed  in  order  to  mount 
the  transducers  and  hydrophones  in  a  3-dimensional 
configuration.   Figure  1  is  a  photograph  of  a  model  of  the 
frame,  showing  the  orientation  of  the  X,  Y,  and  Z  directions 
Figure  2  is  a  schematic  drawing  of  the  sound  fields. 
Figure  3  is  a  photograph  of  the  actual  frame  in  position 
at  the  experiment  site.   The  frame  was  constructed  of 
aluminum  pipe  with  steel  joints.   The  frame  was  designed 
to  be  disassembled  for  portability.   The  fully  assembled 
frame,  less  instrumentation,  weighed  approximately  300  lb. 

b.  Reflective  Noise  Considerations 

Since  it  was  the  objective  of  the  experiment 
to  measure  the  phase  and  amplitude  fluctuation  of  a  sound 
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PHOTOGRAPH  OF  MODEL  OF  FRAME 
FIGURE   1 
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SOURCE 
A:F-33 
B  :  F  -  4 1 
C:F-41 


RCVR 
a:LC-10 
b  :  LC-1 0 

c:LC-10 


PATH  LENGTH 
1 .35  m. 
1 .50  m. 
1 .65  m. 


SCHEMATIC  DRAWING  OF  FRAME,  SHOWING 
LOCATION  OF  SOUND  AXES 

FIGURE   2 
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beam  traveling  a  direct  path  between  the  transducer  and 
the  hydrophone,  it  was  necessary  to  eliminate  any  re- 
flected sound.   Pulse-echo  testing  of  the  frame  configura- 
tion revealed  points  of  specular  reflection  on  the  frame; 
acoustic  absorbant  rubber  was  affixed  to  the  frame  at  these 
points,  and  all  such  reflections  were  reduced  to  at  least 
30  dB  below  the  direct  signal  level. 

c.   Wire  Tension  Considerations 

Scattering  interference  considerations  neces- 
sitated mounting  the  hydrophones  on  small  diameter  wire. 
Such  mounting  caused  vibration  problems,  however,  which 
arose  from  Karman-vortex  shedding  as  water  flowed  past 
the  wires.   Vibration  of  the  mounting  wire  could  cause  an 
accelerometer  effect  in  the  hydrophone  response,  and  more 
importantly,  it  could  cause  a  false  fluctuation  in  the 
received  signal.   This  problem  was  decreased  by  maximizing 
the  tension  of  the  wires. 

2 .   Selection  of  Transducers  (Sources) 

Since  the  experiment  was  to  be  performed  in  a 
frequency  range  from  65kHZ  to  105  kHZ,  a  transducer  with 
a  well  documented  transmit-response  in  this  range,  and  a 
reasonable  degree  of  directivity  was  needed.   One  USRD 
type  F-33  transducer  and  two  USRD  type  F-41  transducers 
were  chosen.   Transmit  responses  for  these  transducers  are: 
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TRANSDUCER         65  kHZ      105  kHZ 
F-33  37dB        48dB 

F-41  37dB        49dB 

dB  reference  1  microbar/volt  at  1  meter 

3.  Selection  of  Hydrophones 

A  hydrophone  with  a  flat  receiving  sensitivity  in 
the  frequency  range  of  interest  was  the  LC-10.   Its 
measured  sensitivity  was: 

6  5  kHZ        105kHZ 
-109dB        -HOdB 
dB  reference  1  volt/microbar 

4 .  Circuit  for  Amplitude  Fluctuation  Detection 
Figure  4  is  a  block  diagram  of  the  instrumentation 

used  to  measure  amplitude  and  phase  fluctuation.   The  phase 
fluctuation  detection  circuit  is  described  in  section  II 
B-5.   Amplitude  fluctuation  detection  is  described  in  the 
following. 

a.   Wave  Synthesizer 

A  General  Radio  Wave  Synthesizer,  model  116 2-A 
was  used  to  generate  a  constant  frequency,  constant  am- 
plitude signal.   The  output' of  the  synthesizer  was  divided, 
part  of  it  going  to  a  Hewlett/Packard  amplifier,  model 
46 7A  and  then  to  the  transducer,  and  part  going  to  the 
phase  meter  as  a  reference  signal  for  the  phase  fluctuation 
measurement.   (section  5) 


22 


to  other   channels 


LC-10 


ONE  CHANNEL  OF  PHASE/AMPLITUDE 
FLUCTUATION  MEASUREMENT  CIRCUIT 

FIGURE    4 
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b.  Transducers 

The  transducers  (F-33  or  F-41)  transmitted 
an  acoustic  signal  through  the  medium  to  the  LC-10 
hydrophone.   One  LC-10  hydrophone  was  used  in  each  of  the 
three  channels,  X,  Y,  and  Z. 

c.  Pre-amplif ier 

The  LC-10  generated  a  low  level  voltage  which 
was  pre-amplif ied  by  a  NUS  constant  gain  pre-amplif ier . 
One  NUS  pre-amplif ier  was  used  in  each  channel.   Two  models 
of  the  NUS  pre-amp.  were  used,  model  2010  which  had  a 
gain  of  approximately  30  dB,  and  model  20  30  which  had  a 
gain  of  approximately  20  dB . 

d.  Filters 

The  signal  was  filtered  by  means  of  a  Krohn- 
Hite  filter.   One  filter  was  used  in  each  channel;  two 
model  3322 's  and  one  model  3342  were  used.   The  signal  was 
band-pass  filtered  at  +300  HZ  on  either  side  of  the  carrier 
frequency.   The  filter  has  a  +20  and  +40  dB  gain  option. 

e.  Envelope  Detector 

The  filtered  signal  was  then  passed  through 
an  envelope  detector,  which  is  shown  schematically  in 
figure  5.   The  envelope  detector  provided  full  wave 
rectification  and  low-pass  filtering  to  remove  the  carrier 
signal.   The  response  of  this  envelope  detector  was  non- 
linear at  low  input  voltages  levels.   That  is,  the  output 
of  the  envelope  detector  was  not  proportional  to  the 
envelope  of  the  modulated  input  voltage  at  low  input 
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voltage  levels.   With  the  aid  of  a  computer,  a  5th  order 
polynomial  curve  was  fitted  to  the  empirically  determined 
response  of  the  envelope  detector.   This  curve  related 
input  AC  voltage,  measured  in  volts  peak-to-peak,  to  out- 
put DC  volts.   One  envelope  detector  was  used  in  each 
direction,  and  a  response  curve  was  computed  for  each, 
f.   Differential  Amplifier 

The  output  of  the  envelope  detector  was  a 
fluctuating  voltage  representing  the  fluctuations  of  the 
incoming  acoustic  signal.   This  envelope  had  a  DC  com- 
ponent which  had  to  be  removed  prior  to  amplification. 
This  DC  voltage  removal  was  accomplished  with  a  dif- 
ferential amplifier,  which  took  the  difference  of  two  in- 
put voltages  and  amplified  it.   One  input  to  the  dif- 
ferential amplifier  was  a  DC  voltage  equal  to  the  DC 
voltage  to  be  removed  from  the  envelope.   The  other  input 
was  the  envelope  itself.   One  differential  amplifier  was 
used  in  each  channel.   Two  of  the  differential  amplifiers 
were  Princeton  Applied  Research  Amplifiers,  model  TM-113, 
which  has  a  variable  output  gain  control,  and  a  band-pass 
filter  option.   Maximum  input  voltage  to  the  PAR  amplifier 
is  0.5  volts.   The  other  differential  amplifier  was  a 
Hewlett/Packard  model  2470A,  which  has  a  variable  output 
gain  control.   Maximum  input  voltage  to  the  Hewlett/Packard 
amplifier  is  20  volts.   The  constant  DC  input  voltage  was 
supplied  by  a  Hewlett/Packard  model  721A  power  supply.   A 
"Heli-pot"  potentiometer  was  used  to  achieve  accurate  and 
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stable  DC  voltages.   A  X10  gain  was  applied  to  the  output 
of  all  three  channels  by  the  variable  gain  of  the  dif- 
ferential amplifiers. 

g.   Tape  Recorder 

The  amplified  envelope  was  then  recorded  on  a 
Precision  Instrument  model  PI-6200  8-track  magnetic  tape 
recorder. 

Figure  6  illustrates  the  progress  of  an  in- 
coming signal  from  the  wave  synthesizer  to  the  tape 
recorder. 

5 .   Circuit  for  Phase  Fluctuation  Detection 
a.   Phase  Meter 

After  the  signal  was  band-pass  filtered,  it 
was  divided,  and  part  of  the  signal  was  fed  into  a  phase- 
meter.   One  phase-meter  was  used  in  each  channel.   Three 
different  types  of  phase  meter  were  used:   Dranetz  Phase- 
Meter  Model  305-PA-3001,  Dranetz  Phase-Meter  Model  305-PA- 
3002,  and  a  Wiltron  Phase-Meter.   The  phase-meter  detected 
the  phase  difference  between  an  input  signal  and  a  reference 
signal  of  the  same  frequency.   The  reference  signal  used 
was  the  input  signal  from  the  wave  synthesizer.   The  out- 
put voltage  of  the  phase-meter  was  proportional  to  the  phase 
difference  between  the  input  and  the  reference  signal.   One 
degree  of  phase  difference  equals  ten  millivolts.   The 
phase  meter's  output  voltage  was  recorded  on  the  8-track 
tape  recorder. 
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C.  LOCATION  OF  EXPERIMENT 

The  experiment  was  carried  out  at  the  Naval  Undersea 
Research  and  Development  Center  Oceanographic  Research 
Tower,  located  approximately  one  mile  off  Mission  Beach, 
near  San  Diego,  California.   The  NUC  Tower  provided  a 
ridigly  fixed  platform  in  the  ocean  which  allowed  con- 
tinuous data  accumulation  and  monitoring.   Figure  7  is  a 
photograph  of  the  tower.   It  is  fixed  by  supporting  pins 
driven  63  feet  into  the  ocean  floor.   Electrical  power 
is  supplied  from  the  shore,  thus  insuring  a  power  supply 
of  relatively  stable  voltages  and  frequency.   The  tower 
is  located  in  approximately  60  feet  of  water  over  a  sand 
bottom.   The  instrumentation  frame  shown  in  Figure  2  was 
mounted  on  a  six  foot  cube  box  standoff  which  was  fixed 
to  a  sled.   Figure  3  is  a  photograph  of  the  frame  mounted 
on  the  sled  on  the  tower.   The  X  axis  was  oriented  in  a 
due  westerly  direction.   The  sled  was  attached  to  a  set 
of  rails  which  extended  from  the  working  area  of  the  tower 
to  the  ocean  floor.   The  frame  was  raised  and  lowered  on 
the  rails  by  means  of  an  electrically  driven  winch,  and 
could  be  positioned  at  any  depth. 

D.  INTERFERENCE  FROM  NEIGHBORING  EXPERIMENT 

A  separate  experiment  was  being  carried  out  at  the  NUC 
tower  in  the  proximity  of  this  experiment.   The  neighboring 
experiment  generated  acoustic  impulses  every  second  from  a 
transducer  positioned  near  the  bottom  beneath  the  tower. 
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These  acoustic  impulses  were  detected  by  the  LC-10  hy- 
drophones of  this  experiment.   The  means  of  reducing  this 
noise  effect  is  described  in  section  IIIb.2. 

E.  ENVIRONMENTAL  CONDITIONS 

The  experiment  was  performed  on  the  8th  and  9th  of 
June,  1972.   The  prevailing  environmental  conditions  are 
listed  in  Table  I.   Of  particular  note  was  the  presence 
of  swell  waves  and  internal  waves.   The  swell  waves  were 
approximately  4  ft,  peak  to  trough,  from  the  southwest, 
and  with  a  period  of  approximately  15  seconds.   The  swell 
waves  persisted  throughout  the  experiment.   The  presence  of 
internal  waves  was  noted  in  the  long  term  temperature 
fluctuations  which  were  observed  at  various  depths  from 
6  June  to  9  June.   These  long  term  temperature  fluctuations 
had  a  period  of  approximately  6-12  minutes.   The  presence 
of  internal  waves  was  also  implied  by  the  varying  depth  of 
the  thermocline  during  the  experiment. 

Tuna  and  other  fish  were  noted  in  the  vicinity  of  the 
tower  during  the  experiment. 

F.  CONDUCT  OF  THE  EXPERIMENT 

1 .   Choice  of  Record  Length 

A  "run"  consisted  of  radiating  constant  amplitude, 
constant  frequency  sound  in  three  orthogonal  directions 
simultaneously  for  approximately  20  minutes.   The  20  minute 
run  time  interval  was  decided  upon  for  the  following  reasons 
Preliminary  tests  had  shown  that  the  sound  amplitude  and 
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phase  fluctuations  had  a  strong  component  frequency  of 
approximately  .06  HZ  (corresponding  to  a  period  of  about 
15  seconds) .   Also,  evidence  of  internal  waves  was  ob- 
served in  the  preliminary  fluctuation  records  in  the  form 
of  long  term  amplitude  and  phase  fluctuation;  these  long 
term  fluctuations  had  periods  of  approximately  6-12  minutes 
A  run  length  long  enough  to  record  the  internal  wave  effect 
would  have  to  include  several  cycles,  minimally  30-40 
minutes.   However,  experiment  time  constraints  and  computer 
processing  limitations  made  it  necessary  to  limit  runs  to 
20  minutes.   This  record  length  would  allow  resolution  of 
the  15  second  surface  wave  effect,  and  identification  of 
the  long  period  internal  wave  effect.   However,  little 
analytical  information  could  be  derived  concerning  the 
latter. 

2 .   Choice  of  Frequencies 

Two  Frequencies  were  chosen  for  this  experiment, 
65  kHZ,  and  105  kHZ.   It  has  been  found  by  Medwin  (1970) 
that  the  bubble  population  in  the  upper  ocean  contains  a 
large  number  of  bubbles  that  are  resonant  at  frequencies 
around  60  kHZ.   In  previous  experiments  by  Rautmann  (1971), 
it  was  found  that  a  peak  dispersion  of  the  sound  speed 
occurs  around  6  5  kHZ;  this  dispersion  was  attributed  to 
the  presence  (and  resonance)  of  bubbles  in  the  medium. 
This  dispersive  effect  results  in  speeds  that  may  be  as 
much  as  8  or  10  meters/second  above  or  below  the  empiri- 
cally predicted  or  velocimeter  measured  values.   Near  this 
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frequency  the  variance  of  the  phase  fluctuation  is  also 
very  large.   At  105  kHZ,  this  dispersive  effect  and  the 
large  variance  of  the  phase  fluctuation  is  negligible. 
It  was  decided  to  run  the  experiment  at  6  5  kHZ 
and  105  kHZ  in  order  to  observe  amplitude  and  phase 
fluctuations  in  a  bubble  affected  and  a  non-bubble  affected 
region. 

3 .  Bathythermograph 

A  bathythermograph  was  taken  during  each  run  to 
obtain  a  representative  gross  temperature  structure  of 
the  medium.   The  BT  drop  was  made  from  the  tower,  ap- 
proximately 40  feet  from  the  frame. 

4 .  Summary  of  Runs 

Runs  were  performed  at  several  depths  between  the 
surface  and  the  bottom,  at  the  two  frequencies,  6  5  kHZ 
and  10  5  kHZ.   These  runs  were  performed  during  the  afternoon 
and  evening  of  8  June  and  the  morning  of  9  June. 

Table  II  is  a  summary  of  the  runs. 
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III.   SIGNAL  PROCESSING 

The  data  gathered  during  the  experiment  were  in  the 
form  of  varying  analog  voltage  records  recorded  on  six 
channels  of  1/4  inch  magnetic  tape  (X,Y,  and  Z  direction 
phase  fluctuations,  and  X,  Y,  and  Z  direction  amplitude 
fluctuations) . 

A.  DIGITIZATION 

The  data  were  digitized  using  a  COMCOR  5000  -  SDS  - 
9300  hybrid  computer  system  of  the  Electrical  Engineering 
Department,  Naval  Postgraduate  School.   It  was  observed 
from  the  time  records  that  the  data  were  generally  band- 
limited  to  less  than  4  HZ.   A  sampling  rate  of  50  HZ 
was  used  in  the  digitization  process  to  minimize  aliasing. 
At  no  time  were  the  fluctuation  records  low-pass  filtered. 
Such  filtering  was  avoided  to  prevent  non-linear  phase 
shifting  of  the  signal's  component  frequencies,  which  could 
result  in  errors  in  the  cross-spectral  analysis.   A  sampling 
pulse  width  of  0.1  millisecond  was  used,  and  the  digitized 
data  were  blocked  in  records  of  512  samples  each.   The 
digitization  process  produced  a  7-track  magnetic  tape  in 
octal  motation,  with  the  six  channels  multiplexed. 

B.  CONVERSION  OF  DATA 

In  order  to  process  the  digitized  data  on  the  Naval 
Postgraduate  School's  IBM  360  computer,  the  7-track  octal 


notation  tape  was  converted  to  9 -track,  hexadecimal  format, 
demultiplexed  tape.   In  this  process  of  conversion,  scaling 
and  "clipping"  were  also  accomplished.   The  conversion 
program  is  listed  on  page  234. 
1.   Scaling  of  Data 

During  the  conversion  process,  it  was  necessary 
to  scale  all  of  the  channels.   The  X,  Y,  and  Z  direction 
phase  fluctuation  channels  were  scaled  by  0.1  to  compensate 
for  the  gain  of  the  digitization  process  and  to  convert 
from  volts  to  degrees  of  phase  difference.   The  scaling 
applied  to  the  X,  Y,  and  Z  direction  amplitude  fluctuations 
was  more  involved.   The  aim  in  this  case  was  to  convert 
the  voltage  records  to  acoustic  pressures  at  the  hydro- 
phones for  each  channel.   This  was  accomplished  by  com- 
puting a  channel  gain  for  each  channel,  based  on  the  gain 
of  each  of  the  equipments.   The  polynomial  curves  described 
in  Section  IIa4.e,  were  employed  here  to  describe  the  response 
of  the  envelope  detectors.   A  channel  was  scaled  as  follows: 
for  a  particular  channel,  the  DC  nulling  voltage  (which 
was  removed  by  the  differential  amplifier)  was  added  to 
the  signal;  the  channel  gain  is  then  applied.   It  was  also 
disired  to  normalize  all  three  channels  of  amplitude 
fluctuation  to  a  common  source  level  and  a  common  trans- 
ducer separation  distance  so  that  comparison  of  channels 
would  be  more  meaningful.   The  above  scaling  -  channel  gain, 
source  level  normalization,  separation  distance  normaliza- 
tion -  was  combined  into  one  overall  gain  factor  which  was 
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applied  to  the  X,  Y,  and  Z  amplitude  fluctuation  channels. 
This  gain  factor  was  different  for  each  channel  and  for 
each  run. 

When  the  analysis  was  completed,  it  became  evident 
that  the  scaling  was  in  error,  and  unfortunately  the 
actual  values  of  the  pressure  levels  at  the  hydrophones 
were  unrecoverable  from  the  data. 
2 .   Clipping  of  Spikes 

In  many  of  the  amplitude  fluctuation  records  there 
were  numerous  voltage  "spikes1,  caused  by  radiation  from 
the  research  project  mentioned  in  Section  IID  which  was 
being  conducted  in  the  proximity  of  this  experiment.   In 
addition  to  this  source  of  noise,  there  were  extraneous 
voltage  spikes  which  either  originated  in  the  electronics 
of  the  equipment,  or  resulted  from  a  physical  impulse  shock 
to  the  system.   Many  of  these  spikes  were  removed  from  the 
time  records  during  the  digitization  process.   Many  of 
those  which  remained  in  the  digitized  time  records  were 
removed  by  means  of  a  simple  "clipper"  routine  in  the 
CONVERT  program.   When  an  extraneous  spike  was  encountered 
in  the  data,  the  previous  non-spike  value  was  substituted 
for  it.   A  threshold  spike  detection  level  had  to  be  es- 
tablished for  the  routine;  the  level  had  to  be  low  enough 
to  detect  the  presence  of  a  spike,  but  large  enough  so 
that  actual  data  would  not  be  interpreted  as  spikes.   It 
was  therefore  impossible  to  remove  all  of  the  spikes  from' 
the  time  records.   Samples  of  the  time  signals  showing 
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the  "raw"  data  and  the  "clipped"  data  are  shown  in  Figures 
8  and  9.   (Note  different  ordinate  scales) 

C.   ANALYSIS  OF  DATA 

1 .   Analysis  Scheme 

Spectral  and  cross-spectral  analysis  of  the  fluc- 
tuation records  was  carried  out  as  follows: 

a.  Mean  Removal  and  De-trending 

The  first  operations  performed  on  each  time 
record  were  removal  of  the  mean  value  of  the  record  and 
the  removal  of  any  linear  trend  of  the  record.   The  mean 
value  (DC  value)  was  removed  so  that  in  the  succeeding 
analysis,  the  fluctuations  of  the  signal  alone  would 
determine  the  nature  of  the  results,  and  so  that  comparison 
of  two  different  records  would  be  meaningful.   Removal 
of  the  linear  trend  of  a  record  was  accomplished  by  sub- 
tracting from  the  time  record  the  least-mean-square  fit 
of  the  record.   This  was  done  to  remove  any  long  term  linear 
change  of  the  signal.   Such  a  change  would  have  resulted 
not  from  the  effects  of  the  medium,  but  from  instabilities 
in  the  electronics  of  the  detecting  circuits.   The  mean 
removal  and  the  de-trending  were  accomplished  by  subroutine 
AVERAGE,  and  TREND  in  the  analysis  package. 

b.  Autocovariance  Function 

An  autocovariance  function  was  computed  from 
a  time  record  by  a  "multiply  and  add"  process  described 
by  : 
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Rx(t)  = Z   X±   Xi+| 

x       N-f   i=l  x        x  T 

where  N  =  total  number  of  points  in  record 
f  =  time  lag 

c.  Power  Sprectral  Density 

Using  the  Wiener-Khinchine  relation,  the  auto- 
covariance  function  was  Fourier  transformed  to  produce  a 
power  spectral  density 

L 

G  (f)  =  4  \   R  (  +  )  cos  (27Tft)  d  + 
x        Jq    x 

where  G  (f)  =  Power  Spectral  Density 
f  -    frequency 
t  =  time  lag 
R(t)  =  Autocovariance  Function 

L  =  Total  time  lag  of  Autocovariance 
A  Parzen  window  was  applied  to  the  auto- 
covariance function  prior  to  transformation  for  smoothing 
of  the  power  spectral  densities. 

d.  The  Cross-Covariance  Function 

A  cross-covariance  function  of  two  channels  of 
interest  (e.g.  phase  and  amplitude  fluctuations  in  the 
same  direction  for  the  same  run,  or  phase  fluctuation  in 
two  different  directions  for  the  same  run)  was  computed. 
Since, R   (+)  =  R   (-+) ,  the  two  parts  of  the  cross- 
covariance  function  were  computed  for  positive  time  delays 
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p   (i)     =   i E   X.  Y. 

xyK    '  N-  +    i=l  l   1++ 


R   (  +  )  =  ™ \        X     Y-;  x-;  n 

yx        N-f    .  ,  1   l+t 

J  i  =  l 

e.   The  Cross-Spectral  Density 

The  cross-spectral  density  was  then  computed 

by  computing  the  cospectral  density  C    and  the  quad- 

xy 

spectral  sensity  Q    as  follows: 


C   (f)  =     2  \   [R   (t)  +  R    (-f)]  cos  (2-rrft)  di 
xy  J   xy       xy 

o 

Q   (f)  =     2\     [R   (f)  -  R    (-  +  )]  sin  (2irft)  di 
xy  J    xy       xy 

o 

A  Parzen  window  was  applied  to  the  cross-covariance  function 
prior  to  transformation  for  smoothing  of  the  resulting 
spectral  densities.   The  cospectral  density  and  the  quad- 
spectral  density  were  combined  to  form  the  cross-spectral 
coherence  function,  and  the  cross-spectral  phase  angle: 


Cx  2 (f)+Q   2 (f) 
Coherence  =    - ^Y. 


Gx(f)    Ky   (f) 


_i  Qxv  (f) 

Phase  Angle  =  Tan    -==* 


C    (f) 
xy 


f.   Display  of  Results 

The  computer  program  used  to  make  the  above 
described  computations  is  listed  on  pages  234  through  252 
Various  functions  are  plotted  by  this  program: 
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Autocorrelation  Function  -  The  autocovariance 
function  normalized  with  respect  to  the  i  =  0 
value 

Cross-Correlation  Function  -  The  cross- 
covariance  function  normalized  with  respect 
to  the  square  root  of  the  product  of  the 
t  =  0  value  of  the  component  autocovariance 
functions 

Auto-Spectrum  Level  -  10Log,0  (power  spectral 
density) ,  in  units  of  dB,  reference  degrees 
for  phase  fluctuation  record.   Relative  dB 
are  used  for  amplitude  fluctuation  records. 
Coherence  Function  -  as  described  in  Section 
d.  above 

Cross-Spectral  Phase  Angle  -  as  described  in 
Section  d.  above 
2 .   Parameters  of  Analysis 
a.   Sampling  Interval 

The  analog  voltage  signals  were  sampled  at 
50  HZ  when  digitized.   This  sampling  rate  was  chosen  to 
insure  that  no  aliasing  of  high  frequencies  occured.   A 
trial  spectral  analysis  of  a  typical  amplitude  fluctuation 
record  revealed  that  the  majority  of  the  energy  of  the 
spectrum  was  at  frequencies  less  than  2  HZ.   In  order  to 
reduce  the  amount  of  computer  core  requirements  and  computer 
time,  the  50  HZ  digitized  data  were  sampled  at  every  8th 
point.   This  process  effectively  sampled  the  original  time 
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signal  at  6.125  HZ,  with  a  sampling  interval  of  0.16 
seconds . 

b.  Resolution 

The  resolution  bandwidth  of  the  power  spectra 
is  related  to  the  autocovariance  lag  time  and  sampling 
interval  by: 

Af    TNLAG) At 

where  NLAG  =  maximum  number  of  time  lags 
At  -    sampling  interval 

c.  Degrees  of  Freedom 

The  number  of  degrees  of  freedom  of  the  power 
spectra  is  defined  as 

n  =  rvr^r^r  /  N  -  sample  size 
NLAG  c 

Degrees  of  freedom  are  an  indication  of  the  normalized 
standard  error  in  the  variance 

error  =  Jl 
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IV.   ANALYSIS  OF  DATA 

A.   RESOLUTION  BANDWIDTH  AND  DEGREES  OF  FREEDOM 

In  computing  each  of  the  autocovariance  functions, 
256  time  lags  were  taken,  each  lag  increment  being  .16 
seconds.   Total  lag  time  was  therefore  equal  to  256X 
.16  =  81.96  seconds.   Therefore,  the  resolution  bandwidth 
of  each  power  spectrum  was 

=  .0122  HZ 


(NLAG)  (At)     (512)  (.16) 

The  sample  size  of  the  time  record  for  the  majority 
of  runs  was  6656  points,  At  =  0.16  seconds.   This  represent- 
ed approximately  17  minutes  of  record.   The  number  of 
degrees  of  freedom  of  runs  of  this  length  was 

2N    _  2(6656)    ?fi 
n    NLAG      512 

and  the  normalized  standard  error  for  the  variance  was 
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=  f^~=    -278 

Runs  that  were  shorter  than  6656  samples  had  correspondingly 
smaller  degrees  of  freedom,  and  are  so  noted  in  the  Summary 
of  Runs,  Table  II. 

B.   OVERVIEW  OF  ANALYSIS  RESULTS 

Six  channels  of  data  (X,  Y,  and  Z  phase  fluctuations, 
and  X,  Y,  and  Z  amplitude  fluctuations)  were  recorded  for 
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5  depths  at  frequency  65  kHZ,  and  4  depths  at  105  kHZ. 
For  various  reasons,  the  data  of  several  of  the  channels 
were  of  poor  quality  or  were  non-existant  when  played 
back.   These  bad  channels  were  due  to  the  presence  of  high 
background  noise,  and  equipment/operator  malfunction. 
The  interpretations  set  forth  in  the  following  sections 
are  based  on  fluctuation  records-  which  were  relatively 
free  of  noise,  and,  in  the  experimenter's  best  judgment, 
reflected  the  effect  of  the  medium  on  the  sound  beam. 

One  can  see  from  Table  II  that  a  large  quantity  of 
data  were  gathered  during  this  experiment.   Time  limit- 
ations prevented  examination  of  many  of  the  different 
aspects  and  implications  of  this  data.   Some  of  the 
more  obvious  features  of  the  data  are  discussed  in  the 
following  sections.   This  coverage  is  by  no  means  extensive, 
but  it  is  adequate  to  point  out  that  much  information 
about  the  medium  can  be  derived  from  a  statistical  study 
of  sound  amplitude  and  phase  fluctuations. 

C.   INITIAL  OBSERVATIONS  AND  INTERPRETATION  OF  SPECTRAL 

ANALYSIS 

The  time  records  and  the  results  of  analysis  of  a 
typical  pair  of  phase  fluctuation  channels  are  presented 
as  figures  10  -  14.   The  two  channels  analyzed  are  the 
phase  fluctuations  in  the  Y  and  Z  direction  of  run  PH-4 
at  a  depth  of  18.8  feet,  with  a  transmission  frequency 
of  65  kHZ. 
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FIGURE   10 
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1 .   Autocorrelation  Function 
a.   Frequency  Content 

In  Figure  11,  observe  that  there  are  three 
predominant  periodicities  in  these  autocorrelation  functions 

(1)  Low  Frequency 

A  very  low  frequency  periodicity  can  be 
observed:  only  about  one  quarter  of  a  cycle  is  included 
in  the  81  second  lag  time.   The  period  of  this  long  term 
fluctuation  which  is  approximately  6  minutes  is  substan- 
tiated by  Figure  13.   It  is  believed  that  this  long  term 
fluctuation  is  the  response  of  the  sound  phase  to  temp- 
erature changes  due  to  internal  waves.   One  can  also  ob- 
serve this  internal  wave  effect  in  the  time  records.   As 
the  medium's  mean  temperature  fluctuates  in  response  to 
the  passage  of  an  internal  wave,  the  sound  speed,  and 
hence  the  phase  also  fluctuates. 

(2)  Surface  Wave  Frequency 

The  most  obvious  periodicity  of  the  auto- 
correlation function  is  the  one  which  recurs  approximately 
every  16  seconds.   This  16  second  periodicity  appears  in 
both  autocorrelation  functions  of  Figure  11,  indeed,  it 
appears  in  all  autocorrelation  functions  of  this  experiment. 
As  stated  before,  the  observed  periods  of  the  surface 
swell  waves  was  also  about  16  seconds.   It  is  believed 
that  orbital  motion  of  the  medium,  in  direct  response  to 
the  surface  wave  action,  caused  temperature  and  salinity 
inhomogeneities  and  bubble  patches  to  pass  into  and  out 
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of  the  sound  field  in  a  periodic  manner.   This  movement 
causes  the  phase  and  amplitude  of  the  sound  beam  to 
fluctuate,  as  explained  in  Section  I.B. 
(3)  High  Frequencies 

There  is  an  unsteady,  low  level  2  to  8  HZ 
fluctuation  superimposed  on  the  autocorrelation  function. 
This  may  be  the  result  of  fine  structure  of  the  medium 
inhomogeneities  or  partly  due  to  random  motion  of  the 
medium,  or  system  noise. 

b.   Initial  Drop  in  Autocorrelation  Function 

There  is  a  sharp  drop  in  the  autocorrelation 
function  at  the  first  or  second  time  lag.   This  initial 
drop  is  due  to  the  presence  of  intermittant  noise  spikes 
in  the  corresponding  time  records  which  are  spaced  more 
than  81  seconds  apart  (the  maximum  time  lag  of  the  cor- 
relation products) .   At  the  +  =  0  value  of  the  autocor- 
relation function  all  spikes  in  the  time  record  are 
multiplied  by  themselves  and  added  to  the  autocorrelation 
sum.   For  any  other  time  lag,  the  spikes  will  not  "line 
up"  again,  and  the  effect  of  the  spikes  will  be  distri- 
buted over  the  entire  autocorrelation  function. 
2  .   Power  Spectrum  Levels 

Figure  12  shows  the  power  spectrum  levels  which 
were  obtained  by  Fourier  transforming  the  autocovariance 
function,  and  then  applying  the  operation  of  10Log,Q  to 
the  transform.   The  6  minute  period  (.003  HZ)  component 
is  less  than  our  resolution  (.01  HZ).   Its  effect  appears 
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but  is  unresolvable  in  the  frequency  range  0  -  0.015  HZ. 
The  strong  frequency  component  due  to  surface  swell  is 
observed  at  approximately  .065  HZ  in  the  spectra  for  both 
the  Y  and  Z  directions.   Components  at  higher  frequencies 
are  at  least  10  dB  down  from  the  spectrum  level  of  the 
f  =  .065  HZ  component,  and  probably  represent  the  effect 
of  other  sea  surface  components,  and  noise  of  the  system. 

3 .  Cross  Correlation  Function 

The  cross-correlation  function  shows  a  strong 
cross-correlation  between  the  two  signals  at  about  16 
second  intervals.   The  envelope  of  the  16  second  period- 
icities suggests  that  a  6  minute  cross-correlation 
periodicity  also  exists.   However,  this  cannot  be  veri- 
fied due  to  the  relatively  short  lag  time. 

4 .  Coherence 

If  the  Y  direction  phase  fluctuation  and  the  Z 
direction  fluctuations  were  completely  unrelated,  the 
coherence  function  would  be  equal  to  zero.   If  the  co- 
herence function  equals  1,  the  two  channels  are  completely 
correlated  and  are  linearly  related.   In  general,  if  the 
coherence  function  is  greater  than  zero  but  less  than 
unity,  one  or  more  of  the  following  three  situations 
exists : 

# 

a.  Extraneous  noise  is  present  in  the  measurements 
of  either  or  both  channels. 

b.  The  relationship  between  the  Z  phase  fluctuation 
and  the  Y  phase  fluctuation  is  not  linear. 
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c.   There  are  inhomogeneities  or  bubble  patches 
affecting  one  channel  which  do  not  effect  the  other  channel; 
that  is,  the  propogation  is  not  isotropic. 

The  coherence  function  shown  in  Figure  14  ranges 
from  0  to  0.75;  no  coherence  function  computed  in  this 
experiment  exceeded  0.8.   Of  the  three  situations  listed 
above,  situation  C  is  the  one  most  likely  to  have  caused 
the  coherence  to  be  less  than  unity. 

As  might  be  expected  from  observing  the  similar- 
ities of  the  autocorrelation  functions,  the  coherence 
function  has  a  relative  maximum  value  at  frequencies  of 
approximately  .005,  .065,  .085,  .125,  and  .170  HZ,  the 
strongest  maximum  being . at  .065  HZ.   These  frequencies 
correspond  to  frequencies  which  are  common  to  the  two 
channels,  the  .065  HZ  frequency  being  the  strongest  common 
frequency. 

5 .   Cross-Spectral  Phase  Angle 

The  cross-spectral  phase  angle  function  has  signi- 
ficance only  at  frequencies  for  which  the  coherence  is 
greater  than  zero.   For  example,  at  .065  HZ,  where  the 
coherence  equals  .75,  the  phase  angle  equals  -30°.   This 
measurement  can  be  explained  in  the  following  way:   Ob- 
serve the  alignment  of  the  16  second  peaks  of  the  two  auto- 
correlation functions,  Figure  11.   Careful  measurement  of 
the  time  lag  value  associated  with  corresponding  peaks 
shows  that  the  Z  peak  always  precedes  the  Y  peak  by  1  to 
1.5  seconds.   In  other  words,  the  correlation  peaks  of  the 
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Z  phase  fluctuation  channel  leads  the  correlation  peaks 
of  the  Y  fluctuation  channel  by  1  to  1.5  seconds.   Es- 
sentially the  same  changes  of  medium  took  place  in  the  Y 
and  Z  directions  but  they  are  happening  1  to  1.5  seconds 
sooner  in  the  Z  direction.   The  frequency  of  this  fluc- 
tuation is  about  .065  HZ;  at  this  frequency,  1.5  seconds 
corresponds  to  approximately  30°.   Therefore  the  cross- 
spectral  phase  angle  is  a  measure  of  the  phase  difference 
between  the  two  channels  in  question,  at  a  particular 
fluctuation  frequency. 

The  fact  that  there  is  a  phase  difference  in  the 
occurance  of  the  phase  fluctuation  of  the  sound  beam  in 
the  Y  and  Z  directions  as  shown  by  the  cross-spectral  phase 
angle,  or  by  direct  measurement  of  the  time  lag  between 
peaks  of  the  autocorrelation  functions,  indicates  that 
moving  inhomogeneities  are  causing  the  16  second  period 
fluctuations.   It  is  suggested  that  it  takes  a  particular 
patch  of  water  on  the  order  of  1  or  1.5  seconds  to  move 
from  the  axis  of  the  Y  direction  sound  field  into  the  axis 
of  the  Z  direction  sound  field,  a  distance  of  approximately 
one  meter;  this  implies  a  speed  of  approximately  65  to 
100  centimeters  per  second.   Fringe  field  interaction 
would  occur  at  about  one-half  meter  at  half  this  speed. 
From  the  swell  wave  point  of  view,  if  the  motion  were 
solely  orbital  at  frequency  0.065  HZ  and  amplitude  of 
about  50  cm.  (this  was  the  amplitude  of  the  observed 


59 


swell) ,  the  orbital  velocity  would  be  approximately  25 
cm/second.  The  moving  patch  explanation  appears  to  be 
quite  reasonable. 

It  was  found  that  a  phase  difference  in  the  oc- 
curance  of  the  16  second  sound  phase  fluctuation  between 
the  Y  and  Z  directions  was  common.   Table  III  shows  the 
cross-spectral  phase  angle  of  the  fluctuations  for  several 
runs.   There  is  consistancy  in  the  magnitude  of  the  angle 
but  not  in  the  sign.   This  can  be  explained  by  the  random- 
ness of  whether,  for  example,  a  hot  patch  starts  from  the 
Y  axis  field  or  the  Z  axis  field  position  at  the  beginning 
of  the  run;  this  position  is  of  course  random,  so  that  the 
sign  of  the  phase  would  be  random. 

D.   CORRELATION  TIME  AS  A  FUNCTION  OF  THE  SOUND  FIELD'S 

RELATION  TO  THE  THERMOCLINE 

For  the  purpose  of  this  work,  correlation  time  is 
defined  as  that  time  required  for  the  envelope  of  the 
autocorrelation  function  to  decay  to  1/e  of  its  initial 
value.   In  cases  where  there  is  an  initial  drop  D  at  the 
first  lag  value,  t  =  .16  seconds,  the  correlation  time  is 
taken  to  be  that  time  required  for  the  envelope  of  the 
autocorrelation  function  to  decay  to  (l-D)/e.   (See 
Section  IV.C.l.b.  for  the  explanation  of  the  cause  of 
steep  initial  drops  in  the  autocorrelation  function.) 

It  is  assumed  that  the  speed  of  sound  in  the  upper 
ocean  is  primarily  a  function  of  temperature  and  bubbles 
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TABLE  III 

Cross  Spectral  Phase  Angle  of  Y  and  Z  direction  Phase 
Fluctuations  at  Fluctuation  Frequency  0.06  HZ 


CROSS-SPECTRAL 
RUN DEPTH     FREO COT PHASE  ANGLE 

FT.      kHZ 


PH-6 

11.3 

105 

.75 

+  30° 

PH-4 

18.8 

65 

.70 

-35° 

PH-3 

25 

65 

.70 

-25° 

PH-8 

26.8 

105 

.55 

-30° 

PH-2 

31 

65 

.25 

+  25° 
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(Rautmann  1971) .   The  periodic  recurrance  of  a  particular 
integrated  speed  of  sound  along  a  sound  axis  implies  a 
stable  (not  changing  with  time)  temperature  microstructure. 
The  correlation  time  of  sound  phase  fluctuations  is  an 
indication  of  the  periodic  recurrance  of  the  integrated 
speed  of  sound  of  the  medium  along  the  sound  axis.   A 
long  correlation  time  indicates  stability  in  the  temperature 
structure  of  the  medium  in  the  vicinity  of  the  sound  field. 
In  other  words,  a  volume  of  water  having  a  slowly  changing 
(stable)  temperature  microstructure  was  moving  through  the 
sound  field  with  periodic  motion  in  response  to  surface 
wave  action  (orbital  motion) .   Short  correlation  time  in- 
dicates instability  in  the  temperature  microstructure; 
that  is,  the  volume  of  water  moving  through  the  sound 
field  periodically  in  response  to  surface  wave  action  has 
a  more  rapidly  changing  temperature  microstructure,  re- 
sulting in  a  less  periodic  recurrance  of  an  integrated 
speed  of  sound  along  the  sound  axis. 

The  behavior  of  the  correlation  times  of  phase  fluc- 
tuations in  the  Y  and  Z  directions  have  been  examined. 
It  is  shown  that  the  correlation  times  of  the  phase 
fluctuations  in  the  Y  and  Z  directions  were  greater 
(greater  stability  of  the  medium)  when  the  sound  field 
was  in  the  thermocline,  than  when  it  was  in  the  well  mixed 
layer. 

Figures  15  to  18  are  diagrams  of  the  termerature 
structure  of  the  ocean  in  the  vicinity  of  the  tower  during 
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runs  PH-2  (31  ft,  65  kHZ) ,  PH-3  (25  ft,  65  kHZ) ,  PH-4 
(18.8  ft,  65  kHZ),  and  PH-5  (7  ft,  65  kHZ).   The  position 
of  the  Y  and  Z  axes  for  each  run  is  also  indicated  on  the 
figures.   This  temperature  structure  varied  during  each 
run,  as  the  thermocline  moved  up  and  down  with  the  medium 
with  a  periodicity  of  6  to  10  minutes.   The  temperature 
structure  shown  is  probably  fairly  representative  of  the 
temperature  that  prevailed  during  the  entire  run,  but 
is  truly  representative  of  the  actual  temperature  struc- 
ture only  at  the  precise  time  at  which  it  was  taken. 
Since  the  run  durations  were  1  1/2  to  3  times  the  internal 
wave  periodicity,  data  processing  gives  average  results 
over  the  changing  medium.   Figures  19,  20,  21,  and  22, 
are  the  17  minute  time  records  of  the  channels  of  interest 
of  runs  PH-2, 3,4,  and  5.   A  star  on  each  time  record  marks 
the  time  the  BT  was  taken.   During  run  PH-5  (7  ft,  65  kHZ) , 
the  Z  axis  was  located  well  within  the  well  mixed  isothermal 
layer.   During  run  PH-4  (18.8  ft,  65  kHZ)  and  Y  and  Z 
axes  were  located  in  the  vicinity  of  the  boundary  between 
the  surface  layer  and  the  first  thermocline.   During  run 
PH-3  (25  ft,  65  kHZ),  and  PH-2  (31  ft,  65  kHZ)  the  Y  and 
Z  axes  were  located  in  a  thermocline. 

The  autocorrelation  function  of  these  runs  are  shown 
in  Figures  23,  24,  25,  and  26.   The  corrected  correlation 
times  are  plotted  as  a  function  of  depth,  and  as  a  function 
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of  the  distance  from  the  Y  axis  sound  field  to  the  lower 
boundary  of  the  surface  mixed  layer.   (See  Figure  27)   It 
can  be  seen  that  the  correlation  times  of  runs  PH-2 
(31  ft,  65  kHZ)  and  PH-3  (25  ft,  65  kHZ)  where  the  Y  and 
Z  axes  were  located  well  within  the  thermocline,  is  greater 
than  the  correlation  time  of  run  PH-5  (7  ft,  65  kHZ) , 
where  the  Y  and  Z  axes  were  within  the  surface  layer,  and 
run  PK-4  (18.8  ft,  65  kHZ),  where  the  Y  and  Z  axes  were 
located  in  the  transition  region  between  the  surface  layer 
and  the  thermocline.   It  is  believed  that  this  behavior 
of  the  correlation  time  is  due  to  the  following. 

Measuring  the  decay  time  of  the  autocorrelation  function 
envelope  to  (l-D)/e,  vice  the  decay  time  of  the '16  second 
periodicities,  measures  the  de-correlation  of  the  long- 
term  (6  to  10  minute  periodicities)  phase  fluctuation. 
The  long  term  fluctuations  are  believed  to  result  from 
long  term  temperature  fluctuations  of  the  medium  of  the 
sound  field,  which  reflect  the  medium  movement  in  response 
to  the  passage  of  an  internal  wave. 

In  the  surface  mixed  layer,  the  temperature,  salinity 
and  bubble  microstructure  is  random;  Skudrzyk  (1963)  de- 
scribes this  random  distribution  as  "temperature  patches" 
of  varying  sizes.   On  the  other  hand,  since  the  work  of 
Woods  (1968) ,  a  picture  of  the  more  stable  thermocline 
has  evolved.   The  thermocline  is  said  to  consist  of 
alternate  sheets  of  high  temperature  gradients  (order  of 
cm.  thick)  and  plates  of  low  temperature  gradient  water 
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(order  of  a  meter  thick),  both  parallel  to  the  sea  surface, 
If  a  volume  of  this  well  structured  thermocline  water  is 
moved  vertically  into  and  out  of  a  sound  field  with  an 
approximate  periodicity,  as  in  response  to  the  passage  of 
an  internal  wave,  the  sound  speed  fluctuation,  and  there- 
fore the  phase  fluctuations  should  reflect  the  periodic 
recurrence  of  particular  temperatures.   In  contrast  to 
this,  consider  the  surface  layer  where  the  temperature 
structure  is  random.   There  will  be  only  random  recurrence 
of  particular  temperature  structures  for  the  same  long 
term  medium  movements  since  the  motion  of  the  well  mixed 
layer  has  a  superimposed  random  horizontal  and  vertical 
motion  of  the  inhomogeneities .   Hence,  the  long  term 
(envelope)  fluctuations  of  the  autocorrelation  function 
will  decorrelate  faster  (will  have  a  shorter  correlation 
time)  in  the  surface  layer  than  in  the  thermocline.   This 
behavior  is  in  fact  observed  in  the  shorter  correlation 
time  of  runs  PH-5  which  was  located  in  the  surface  layer, 
and  PH-4  which  was  located  in  the  transition  region  be- 
tween the  surface  layer  and  the  first  thermocline,  and  the 
longer  correlation  times  of  runs  PH-3  and  PH-2,  both  of 
which  are  located  in  a  thermocline. 

E.   Z  DIRECTION  CORRELATION  TIME  VERSUS  Y  DIRECTION  CORj 

RELATION  TIME 

In  runs  PH-2  (31  ft,  65  kHZ)  and  PH-3  (25  ft,  65  kHZ) 
which  took  place  in  the  thermocline,  the  Z  direction  phase 
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fluctuations  have  a  longer  correlation  time  than  the  Y 
direction  phase  fluctuation;  while  in  run  PH-4  (18.8  ft, 
65  kHZ)  which  took  place  in  the  mixed  layer,  the  Y  and  Z 
direction  phase  fluctuations  have  approximately  equal 
correlation  times.   This  may  well  be  due  to  the  presence 
of  stable  layered  patches  of  isothermal  water  in  the  thermo- 
cline, and  isotropically  random  patches  of  isothermal 
water  in  the  mixed  layer.   As  stated  in  Section  IV. D,  by 
defining  the  correlation  time  to  be  the  time  associated 
with  the  envelope  decay,  a  correlation  time  is  a  measure 
of  the  stability  of  the  medium  over  long  term  changes; 
in  this  case  changes  resulting  from  internal  waves.   Long 
term  periodicities  of  a  phase  fluctuation  depend  on  a  par- 
ticular value  of  integrated  sound  speed  recuring  in  the 
sound  field  at  regular  (long  term)  intervals.   As  previously 
mentioned,  there  is  evidence  that  there  are  regions  of 
thermocline  "sheets"  and  "plates"  parallel  to  the  surface 
of  isothermal  water  in  the  thermocline. 

If  the  medium  undergoes  vertical  movement,  the  extent 
of  the  periodicity  of  long  term  phase  fluctuations  in  the 
Y  direction  will  depend  on  a  recurring  isothermal  plate 
at  the  Y  axis,  since  the  Y  axis  sound  beam  travels  hori- 
zontally through  the  isothermal  plate.   In  contrast,  the 
Z  axis  sound  beam  travels  vertically  through  horizontal 
plates  of  isothermal  water.   The  integrated  sound  speed  over 
the  path  is  the  average  temperature  of  the  horizontal 
plates.   Long  term  periodicities  of  the  phase  fluctuation 


79 


in  the  Z  direction  depend  on  the  recurrence  of  an  average 
sound  speed  over  one  or  more  isothermal  plates;  such  a 
periodicity  in  the  Y  direction  depends  on  the  recurrence 
of  a  particular  isothermal  plate. 

Therefore,  in  a  thermocline,  where  a  fairly  stable 
gradient  prevails,  long  term  phase  fluctuations  are  more 
stable  in  the  Z  direction  than  in  the  Y  direction  because 
the  orbital  movement  of  the  medium  is  less  likely  to  cause 
deviations  from  a  recurring  average  sound  speed  in  the  Z 
direction  than  in  the  Y  direction. 

Note  that  for  run  PH-4  (18  ft,  6  5  kHZ)  which  was  in 
the  transition  region  between  the  thermocline  and  the  mixed 
layer,  the  Y  and  Z  correlation  times  are  nearly  equal. 
The  mixed  layer  is  composed  of  randomly  distributed  patches 
It  is  therefore  equally  likely  to  have  an  integrated  sound 
speed  repeated  in  the  Y  direction  as  in  the  Z  direction. 

Thus,  for  the  cases  examined,  sound  phase  correlation 
times  are  not  isotropic  in  the  thermocline,  but  are  nearly 
isotropic  in  the  mixed  layer. 

F.   Y  DIRECTION  VARIANCE  VERSUS  Z  DIRECTION  VARIANCE 

Figure  2  8  is  a  graph  of  the  variance  of  Y  and  Z 
direction  phase  fluctuations  for  runs  PH-2  (31  ft,  65  kHZ) , 
PH-3  (25  ft,  65  kHZ),  and  PH-4  (18.8  ft,  65  kHZ).   The 
variance  of  a  time  record  is  defined  as  the  zero  lag  value 
of  the  autocovariance  function.   Because  of  the  sharp 
initial  drop  in  many  of  the  autocovariance  functions  due 
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to  noise  spikes,  the  value  of  the  autocovariance  function 
at  +  =  0.16  or  0.32  seconds  is  a  more  accurate  estimate 
of  the  variance  than  the  f  =  0  value.   Therefore,  in 
comparing  the  variances  of  the  Y  and  Z  direction  fluctua- 
tions of  runs  PH-2,3,  and  4,  this  corrected  variance  was 
used. 

Chernov  (1960)  has  shown  that  the  variance  of  a  phase 
or  amplitude  fluctuation  in  a  random  medium  is  proportional 
to  the  length  of  the  transmission  path.   The  Y  and  Z  dire- 
ction paths  were  of  different  lengths,  and  although 
Chernov's  assumptions  are  not  fulfilled  in  the  thermocline, 
all  variances  were  normalized  to  one  meter  by  dividing 
each  variance  by  its  path  length.   The  corrected  and  nor- 
malized variance  versus  run  depth  is  plotted  in  Figure  28. 
The  Y  direction  variance  is  always  greater  than  the  Z 
direction  variance.   This  may  also  be  due  to  the  layered 
structure  of  the  medium.   If  horizontal  plates  of  various 
thicknesses  are  moving  through  the  sound  field  due  to 
orbital  motion,  the  integrated  phase  fluctuations  will  be 
more  sensitive  to  a  vertical  movement  in  the  Y  direction 
than  in  the  Z  direction,  because  the  Y  direction  sound 
beam  passes  horizontally  through  a  horizontal  isothermal 
"plate".   A  small  movement  could  shift  a  new  "plate"  and 
therefore  a  sharp  temperature  change  into  the  Y  axis, 
while  the  integrated  speed  and  temperature  over  the  Z 
direction  would  not  be  so  greatly  effected. 
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The  Z  direction  variance  of  the  phase  fluctuations  for 
run  PH-5,  which  took  place  in  the  surface  mixed  layer,  is 
much  greater  than  the  variances  of  runs  PH-2,  PH-3,  and 
PH-4 .   This  is  an  expected  result  due  to  two  characteristics 
of  the  surface  mixed  layer:  the  random  orientation  of 
temperature  patches  and  the  presence  of  bubble  populations 
which  have  been  found  to  be  resonant  at  65  kHZ. (Rautmann, 
1971,  Fitzgerald,  1972) 

It  is  concluded  that  the  variance  of  the  phase  fluc- 
tuation in  the  thermocline  is  not  isotropic. 

G.   CORRELATION  OF  PHASE  AND  AMPLITUDE  FLUCTUATIONS 

Chernov  (1950)  has  derived  expressions  for  the  temporal 
autocorrelation  function  of  the  phase  and  amplitude  fluc- 
tuations of  a  sound  beam  in  an  isotropic,  stationary  random 
medium.   He  concludes  that  the  temporal  autocorrelation 
function  of  the  phase  fluctuations  equals  the  temporal 
autocorrelation  function  of  the  amplitude  fluctuations. 
The  autocorrelation  function  of  phase  and  amplitude  fluc- 
tuations of  this  experiment  can  be  compared  to  see  how 
closely  they  approach  the  Chernov  equality. 

In  examining  the  autocorrelation  functions,  it  will 
be  noted  that  in  many  cases  the  autocorrelation  function 
has  a  very  steep  initial  drop  at  t  =  .6  seconds  after  which 
it  decays  more  normally.   It  is  believed  that  this  sharp 
initial  drop  is  due  to  extransous ,  widely  spaced  noise 
spikes,  as  explained  in  Section  IV  C.l.b.   Therefore, 
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in  comparing  phase  and  amplitude  fluctuations,  the  effects 
of  the  extraneous  spikes  should  be  removed  from  each  func- 
tion prior  to  comparison.   This  could  be  done  by  taking  the 
0.16  second  or  0.32  second  value  of  the  autocovariance 
function  as  the  normalizing  factor  in  each  case.   This 
re-calculation  of  the  autocovariance  function  was  not  done, 
for  lack  of  time,  but  a  roughly  equivalent  correction  was 
made,  that  of  merely  subtracting  the  initial  drop  from  each 
of  the  autocorrelation  functions.   Figure  29  shows  the 
autocorrelation  function  of  the  Z  direction  phase  and 
amplitude  fluctuations  from  PH-8  for  65  kHZ,  26.8  ft.' 
depth.   Figure  30  shows  the  same  two  autocorrelation  func- 
tions with  the  initial  drop  subtracted  from  each  function. 
It  can  be  seen  that  the  phase  and  amplitude  fluctuation 
autocorrelation  functions  are  very  similar  in  this  case. 
This  is  an  example  of  excellent  agreement  between  the  phase 
and  amplitude  fluctuation  autocorrelation  functions.   Figure 
31  shows  the  X  direction  phase  and  amplitude  autocorrela- 
tion functions  from  PH-7  for  the  65kHZ,  4  4.2  ft.  depth  run. 
In  this  case,  the  phase  fluctuation  is  greatly  affected 
by  extraneous  noise  spikes,  and  by  high  frequency 
periodicities.   The  extraneous  noice  spikes  cause  the  initial 
sharp  drop  at  t  =  .16  or  .32  seconds.   The  high  frequency 
periodicity  is  either  high  freqnency  noise,  or  the  effects 
of  aliasing  of  high  frequency  noise.   In  either  case, 
these  high  frequency  signals  are  not  caused  by  the  effects 
of  the  ocean  medium  on  the  sound  wave,  as  explained  in 
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Section  IV.C.l.b;  hence,  these  high  frequency  effects  are 
not  considered  in  this  analysis.   In  Figure  32 ,  the  ex- 
traneous noise  spike  effect  has  been  removed,  and  the  phase 
autocorrelation  function  has  been  smoothed.   It  can  be  seen 
that  the  phase  and  amplitude  autocorrelation  functions  are 
still  similar  in  appearance,  and  have  relative  maxima  and 
minima  at  about  the  same  lag  times. 

Other  pairs  of  phase  and  amplitude  autocorrelation 
functions  are  shown  in  Figures  33  -  38.   If  the  initial 
sharp  drop  is  ignored,  it  will  be  observed  that  their 
behaviors  are  essentially  the  same.   Thus,  for  the  case 
observed  here,  it  appears  that  the  Chernov  statement 
holds,  to  varying  degrees.   The  reasons  why  the  Chernov 
equality  fails  are:  1)  the  presence  of  different  amounts 
of  extraneous  noise  spikes  in  each  signal  causing  a  sharp 
initial  drop,  and  2)  the  presence  of  extraneous  high 
frequency  signals  in  one  of  both  channels  (noice) ,  and 
most  importantly,  3)  the  phase  fluctuations  are  sensitive 
to  long  term  temperature  changes  while  the  amplitude  fluc- 
tuations are  not.   Chernov  assumed  that  random  inhomogeneities 
in  the  sound  field  were  the  sole  cause  of  fluctuations; 
he  did  not  consider  long  term  anisotropic  temperature 
changes,  such  as  those  that  result  from  the  passage  of  an 
internal  wave. 
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V.   CONCLUSIONS 

The  following  conclusions  are  drawn  based  on  the 
preceeding  analysis: 

A.  Phase  and  amplitude  fluctuations  of  a  65  kHZ  and  105 
kHZ  sound  beam  in  the  surface  layer  of  the  ocean  are  caused 
by  the  motion  of  medium  inhomogeneities  of  varying  acoustic 
impedance  in  the  sound  field.   Inhomogeneity  motion  arises 
from  two  main  sources : 

1.  Orbital  motion  of  the  medium  in  response  to  surface 
wave  action,  and 

2.  Motion  of  the  medium  in  response  to  the  passage  of 
internal  waves. 

The  dominant  frequency  components  of  the  phase  and 
amplitude  fluctuations  are  those  of  the  surface  waves  and 
internal  waves . 

B.  The  envelopes  of  the  autocovariance  functions  of  the  phase 
fluctuations  in  a  horizontal  and  a  vertical  sound  beam  of 

65  kHZ  decorrelate  (reduce  to  1/e  of  their  initial  value) 
faster  in  the  surface  mixed  layer  than  in  the  thermocline. 
It  is  believed  that  this  is  an  indication  of  two  character- 
istics of  the  mixed  layer:   the  greater  instability  of  the 
temperature  structure  in  the  mixed  layer,  versus  that  in  the 
thermocline,  and  the  presence  of  bubbles  resonant  at  65  kHZ 
in  the  mixed  layer. 
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C.  In  the  thermocline,  the  autocovariance  function  of  the 
phase  fluctuations  of  a  horizontal  sound  beam  decorrelates 
faster  than  the  autocovariance  function  of  the  phase  fluc- 
tuations of  a  vertical  sound  beam.  It  is  believed  that 
this  results  from  a  layered  "plate"  temperature  structure 
of  the  thermocline.   Therefore,  the  correlation  time  in 

a  thermocline  is  not  isotropic. 

D.  In  the  thermocline,  the  variance  of  the  phase  fluctua- 
tion of  a  65  kHZ  sound  beam  in  the  horizontal  direction  is 
greater  than  the  variance  of  the  phase  fluctuations  in  the 
vertical  direction.   It  is  believed  that  this  is  a  result 
of  a  layered  "plate"  temperature  structure  of  the  thermo- 
cline.  Therefore,  the  variance  in  the  thermocline  is  not 
isotropic . 

E.  From  C.  and  D.  above,  it  is  concluded  that  the  medium 
in  the  thermocline  is  anisotropic. 

F.  The  autocorrelation  function  of  the  phase  fluctuation 
approximately  equals  the  autocorrelation  function  of  the 
amplitude  fluctuations  for  a  65  kHZ  sound  beam  traveling 
through  1  meter  of  the  surface  layer  of  the  ocean.   This 
result  reinforces  and  extends  the  theoretical  work  of 
Chernov  (19  50)  which  predicted  that  the  autocorrelation 
functions  of  the  phase  and  amplitude  fluctuations  of  a 
sound  beam  traveling  through  a  random,  stationary  isotropic 
medium  would  be  equal. 
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VI.   RECOMMENDATIONS  FOR  FUTURE  EXPERIMENTS 

A.   EXPERIMENTAL  TECHNIQUE 

The  results  of  this  experiment  demonstrate  the  feasi- 
bility of  obtaining  useful  information  about  the  micro- 
structure  of  a  volume  of  ocean  water  by  studying  the  sound 
phase  and  amplitude  fluctuations  of  sound  passing  through 
the  volume.   It  also  points  out  some  of  the  pitfalls  in 
making  such  measurements.   For  those  who  repeat  this 
experiment,  or  perform  similar  experiments  in  the  future, 
the  following  recommendations  are  offered. 

1.  Use  differential  amplifiers  that  are  not  limited 
to  a  low  input  voltage.   This  will  allow  operation  well 
above  many  of  the  noise  levels  of  the  system.   The  PAR  113 
differential  amplifier,  although  ideal  from  the  point  of 
view  of  internal  electrical  noise,  was  poorly  suited  for 
this  purpose. 

2.  To  obtain  any  variance  information  from  the  am- 
plitude fluctuation  records,  the  system  gains  must  be  known 
exactly  (phase  fluctuation  measurements  will  not  be  effected 
by  varying  gains  or  voltage  levels) .   Careful  measurement 

of  these  gains,  preferably  before  the  experiment,  is  the 
only  way  to  obtain  such  information.   System  gain  measure- 
ments should  be  done  with  the  detection  circuit  fully  as- 
sembled, vice  piecemeal. 
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B.   AREAS  FOR  FUTURE  INVESTIGATION 

1.  Much  more  information  can  be  gleaned  from  the  data 
of  this  experiment.   Analysis  results  of  runs  not  pre- 
viously considered  are  included  as  Appendix  A.   It  is  felt 
that  further  analysis  of  these  data  would  be  fruitful. 

2.  Concurrent  measurement  of  sound  amplitude  and  phase 
fluctuation,  particle  velocity,  and  temperature  in  the 
sound  field  would  shed  more  light  on  the  relationship 
between  medium  motion  and  sound  phase  and  amplitude  fluc- 
tuation.  By  making  such  concurrent  measurements  of  particle 
velocity,  temperature,  and  sound  phase  and  amplitude 
fluctuation,  it  should  be  possible  to  deduce  the  role  of 
resonant  bubbles  in  the  acoustic  fluctuations. 
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APPENDIX  A 
Analysis  Results  For  Various  Runs 

Table  II  lists  the  experimental  runs  which  were 
analyzed  by  the  computer  programs  CONVERT  and  ANALYSIS . 
Various  outputs  of  this  analysis  are  displayed  throughout 
the  body  of  the  thesis. 

Appendix  A  contains  the  remainder  of  the  analysis. 
This  information  is  included  to  permit  comparison  with 
future  studies  in  this  field,  and  for  further  analysis 
of  this  data  at  some  future  time.   A  typical  "set"  of 
results  for  a  given  run  are  (1)  Temporal  Autocorrelation 
Function,  (2)  Power  Spectrum  Level,  (3)  Temporal  Cross- 
Correlation  Function,  and  (4)  Coherence  Function  and 
Cross-Spectral  Phase  Angle.   Some  sets  are  not  complete, 
as  portions  of  them  were  displayed  in  the  body  of  the 
thesis . 

Interpretation  of  these  outputs  is  explained  in  Section 

IV-C.   The  ordinate?  of  the  cross-correlation  function 

graphs  are  labeled  incorrectly;  they  should  be  labeled  as 
in  Figure  13. 
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Y-SCRLE--=2.00E+02  UNITS  INCH. 

CROSS   SPECTRAL   PHASE   RNGLE   K-PM,Z-PM 

RUN  PH-3 


128 


K-5CRLF.----2.0QE+Q1  UNITS  INCH. 
Y-SCflLF-2-  OQF-OJ  UNITS  INCH. 

TEMPORAL  AUTOCORRELATION  FN, X-PM, V-PM 
RUN  PH-3,  FILE  6  OF  C0N6 


129 


032 


024 


_Y-PK 


S-SCnLE-i.OOE-Oi  UNITS  INCH. 
Y-SCRLE=1.00E+01  UNITS  INCH. 

PQRiER   SPECTRUM  LEUEL 
RUN   PH-3,    FILE   6   CF   CGN6 


GB1    K-PM,V-PM 


130 


r-q 


031 


X-SCflLE->i.00t>00  UNITS    INCH, 
Y-SCHLf>5 . 00E-02  UNITS    INCH. 

CROSS-CORRELRTIGN  FN,X-PM,Y-PM 

RUN  PH-3,  FILE  6  OF  CQN6 


131 


,r_:Hn 


QS4 


X-SCRLE----1.00E-0J  UNITS  INCH 
V-SCRLE=1.00E+00  UNITS  INCH 

CGHERENCE  FUNCTION 


X-P 


M     U 


PM 


X-SCRLE-l.OQE-Oi   UNITS    INCH. 
Y-SCfiLE^2-00E+02   UNITS    INCH. 

CROSS  SPECTRAL  PHRSE 


RNGLE  K-PM.V-PM 


RUN  PH-3 


132 


yV- 


_7-PM 


JC2 


r7-R« 


X-SCRLE=2.00E+01  UNITS  INCH. 
y-SCRLE--=2.00E-01  UNITS  INCH. 

TEMPORAL   RUTQCQRRELflTION   FN.Z-PM.Z-F 

RUN   PH-3,    FILE   6   OF   CQN6 


133 


f-ni 


n  t  3 
u  u  .1 


j:-pm 


K-SCRLE-l.OOE-01  UNITS. INCH 
V-SCRLE-l.OOE+01  UNITS  I 


PORER   SPECTRUM   LEUEL    EDB)    Z-PM,Z-RM 
RUN   PH-3,    FILE   6   GF   CCN6 


134 


K-SCflLE-l.OOE+OQ  UNITS  INCH. 
Y-SCRLE->5.00E  -02  UNITS  INCH- 


CROSS 
RUN 


PH-3 


CORRELATION  FN,Z-PM 

FILE  6  OF  C0N6 


7 

i  L 


RM 


135 


oc^- 


K-SCRLE=i-OOE-01  '^NITS    I 
Y-SCP,LE==i-OOE+00  UNITS    INCH 

COHERENCE  FUNCTION 


Z-PM,Z-F 


A 


K-SCRLE-l.OOE-Ol  UNITS  INCH. 
Y-5CflLE--=2.00E+02  UNITS  INCH. 

CROSS  SPECTRAL  PHRSE 


RNGLE  Z-PM 


?    L. 


RM 


136 


X-SCRLE-l.OOE-01 
Y-SCRLE=1-00E+01 


rs 


UNI 
UNITS 


POWER  SPECTRUM  LEUll 


INCH. 
INCH. 

r 


RUN   PH-4,    FI 


1  LL 


r 


OF 


DB 

;gnb 


Z-PM,2-RM 


137 


K-SCRLF.-l- 00E->00  UNITS  INCH. 
V-SCflLE-l.OOE-01  UNITS  INCH- 

CROSS-CORRELn.TION  FN,Z-PM,Z-flM 
RUN  PH-4,  FILE  7  OF  C0N6 


138 


DC4 


X3HR 


K-SCRLE=1-00E-01  UNITS  INCH 
Y-SCRLE=l-OOE+00  UNITS  INCH 

COHERENCE  FUNCTION 


Z-PM.Z-RH 


X-SCRLE-l-OOE-Ol  UNITS  INCH. 
V-SCRLE----2.00E+02  UNITS  INCH. 

CROSS  SPECTRRL  PHRSE  ANGLE, Z-PM, Z-RM 

RUN  PH-4 


139 


V-PM 


NAf^/ 


X-SCflLE^--2.00E+0i   UNITS    INCH. 
Y-SCRLE^.OOE-Oi   UNITS    INCH. 


TEMPORRL  RUTOCORRELRTION  FN,X-PM,Z- 
RUN  PH-4,  FILE  7  OF  C0N6 


140 


-PM 


0:2 


03.3 


DIM 


K-SCRLE=i.00E-01'  UNITS  INCH. 
Y-SCRLE=1. OOE+OI  UNITS  INCH. 

POWER   SPECTRUM   LEUEL    (DB) 
RUN   PH-4,    FILE   7   OF   C0N6 


J!-PJ1 


vi  "7  ■  D  M 


X-PM,Z-P 


141 


o 

-I 

C) 


r 


i     tCSR 


032 


03-1 


K-SCHLE----1.00E+00  UNITS    INCH. 
V-SCflL£:---5.00t:-02  UNITS    INCH- 

CROSS-CORRELATION   FN'.X-PM,Z-PM 
RUN   PH-4,    FILE   7   OF   C0N6 


142 


sCk^ 


.Acfl 


0D4 


X-SCRLE-i-OOE-Ol  UNITS  INCH 
Y-SCRLE-i.OOE+00  UNITS  INCH 

COHERENCE  FUNCTION 


X-P 


M 


PM 


PHP.S 


X-SCFILE-1.0QE-01   UNITS    INCH. 
Y-SCRLE^2-00E+02  UNITS    INCH- 


CROSS  SPECTRAL  PHASE  RNGLE, X-PM, Z-PM 

RUN  PH-4 


143 


^PW^\^ 


3 


K-SCRLE-2.00E+01  UNITS  INCH. 
Y-SCFILE----2.00E-Q1  UNITS  INCH. 

TEMPORAL  RUTQCQRRELRTIQN  FNpX-RMpV-Rfl 

RUN  PH-4,  FILE  7  OF  C0N6 


144 


occ 


o:'2 


0C3 


rj:.M 


K-SCRLE=1.00E-01  UNITS  INCH. 
V-5CRLr---2.00t+01  UNITS  INCH. 

PONER  "SPECTRUM   LEUEL    (DB) 
RUN   PH-4,    FILE   7   OF   CGN6 


X-RM,Y-RM 


145 


9    -v, 


032 


004 


jDCOR 


X-SCflLE^LOOE+OO   UNITS    INCH. 
V-SCRLE^5. OOE-02  UNITS    INCH. 

CROSS-CORRELRTION   FN,X-RM,Y-RM 

RUN   PH-4.    FILE   ?   OF   GQN6 


146 


X3HF1 


K-SCRLE=1. OOE-01  UNITS  INCH 
Y-SCRLE-1. OOE+QO  UNITS  INCH 

COHERENCE  FUNCTION 


K-RMpV-RM 


K-SCRLE-1. OOE-01  UNITS  INCH. 
V-SCRLE-2.00E+02  UNITS  INCH. 

CROSS  SPECTRAL  PHRSE 


RNGLE,X-RM,Y-RM 


RUN  PH-4 


147 


JC2 


o:< 


X-SCALE>2--00F+01  UNITS  INCH 
Y-SCRLF--=2.QQE-Q1  UNITS  INCH 
T 


EMPORRL   RUTOCORRELRTION 
PH-4,    FILE   7   OF   C0N6 


tr 


M 


X-RM,Z-RM 


148 


0C2 


P_. 


OOE-01  UNITS  INCH. 
Y-SCRLE=2.00F+01  UNITS  INCH. 

POWER  SPECTRUM  .LEUEL  (DB) 
RUN  PH-4,  FILE  7   Or  CGN6 


X-RhUZ-RM 


149 


X-SCflLE=i . OOE+00  UNITS  INCH. 
V-SCRLE-l.OOE-01  UNITS  INCH. 


CROSS-CORRELRTION  FN, X-RM, Z-RM 
RUN  PH-4,  FILE  7  OF  C0N6 


150 


dm 


-yXSHn 


K-SCRLE--=1.00E-;-01  UNITS  INCH- 
Y-SCRLE-1.  OOE+00  UNITS  INCH. 

COHERENCE   FUNCTION 


X 


RM 


R 


■vH 


X-SCRLE-l.OOF-Oi   UNITS   INCH. 
V-SCRL£>2.0Qt>02   UNITS    INCH. 

CROSS   SPECTRRL   PHASE 


RNGLE,X-RM,Z-RM 


RUM  PH-4 


151 


"^^^^^^ 


K-5CRLE=2-00E+01  UNITS  INCH, 
Y-SCRLE--2.  00E-Q1  UNITS  INCH. 

TEMPQRRL  RUTQCGRRELRTIGN  FN,Y-RM,Z-RM 

RUN  PH-4?  FILE  7  OF  C0N6 


152 


K-SCRLE=1.00E-01  UNITS  INCH. 
V-SCRLE---1. OOE+Oi  UNITS  INCH. 

PGHER  SPECTRUM  LEUEL  ( DB 
RUN  PH-4,  FILE  7  OF  C0M6 


v-RM,z-f 


153 


X-SCRLE-l.QOE+00  UNITS  INCH. 
Y-SCRLE=1.00E-01  UNITS  INCH. 

CROSS-CORRELRTION  FN,V-RM,2-RM 
RUN  PH-4,  FILE  7  OF  CGN6 


154 


K-SCRLE-l-OOE-Ol  UNITS  INCH. 
Y-SCRLE-LOQE+OQ  UNITS  INCH. 

COHERENCE  FUNCTION  V-RM,Z-RM 


K-SCRLE----1-OOE-01  UNITS  INCH, 
V-SCRLE--?2-00E+02  UNITS  INCH. 

CROSS  SPECTRRL  PHRSE  RNGLE, Y-RM,Z-RM 

RUN  PfM 


155 


:-PM 


V^-PK 


K-SCRLE=2-00E+01  UNITS  INCH. 
Y-SCRLE=2.00E-Q1  UNITS  INCH. 

TEMPORRL  RUTOCQRRELRTIQN  FN,X-PM 

RUN  PH-4,  FILE  7  OF  CGN6 


V-PM 


156 


a 


~"\.x~pm 


_v-PM 


X-SCRLEyi.tlOE-01 

y-scnLE==i-ooE+oi 


UNITS  INCH. 
UNITS  INCH. 


POWER  SPECTRUM  LEUEL  (DB)  K-PM.V-PM 
RUN  PH-4,  FILE  7  QF  C0N6 


157 


ICC3P 


QQ< 


K-SCRLE----1.00F  +  00  UNITS  INCH- 
V-SCRLE-5.00E-02  UNITS  INCH. 

CROSS-CORRELnTIQN  FN, X-PM, Y-PM 

RUN  PH-4,  FILE  ?  OF  CGN6 


158 


X-SCRLE---\L.0QE-01   UNITS    INCH 
V-SCRLE=--l-OOE+00   UNITS    INCH 

COHERENCE   FUNCTION 


I   I      f~\   I    I —  I   J  /~V  I   (  i — ^  i — 


K-PMpV-PH 


X-SCRLE=1.00E-01   'JNITS    INCH. 
V-SCRLE-2.00E+02  UNITS    INCH. 

CROSS   SPECTRAL   PHRSE 


RNGLE,X-PM,Y-PM 


RUM    PH-1 


159 


K-SCRLF.-2.0QE+01  UNITS  INCH. 
^-SCRLE-2.00E-01  UNITS  INCH. 

TEMPQRRL   RUTQCORRELRTION   FN.V-P 


k 


V-RM 


RUN   PH-4 


4:  y 


r 


ILE   7   OF   C0N6 


160 


y-m 


y-pe. 


K-SCRLE=1-00E-01  UNITS  INCH. 
V-SCRLE--5.00E+00  UNITS  INCH- 

POUER   SPECTRUM   LEUEL    ( DB1 

RUN   PH-4,    FILE   7   OF   CGN6 


;J-Pt1 


Y 


RH 


J.61 


X-SCRLF>=i.OOE"+00  UNITS  INCH- 
V-SCRLE-S.00E-02  UNITS  INCH. 

CROSS-CORRELATION  FN,Y-PM,Y~nM 
RUN  PH-4,  FILE  7  OF  C0N6 


162 


JtfHn 


K-SCflLE=l.Q0E-01  UNITS  INCH 
V-SCRLE----1-00E+00  UNITS  INCH 

COHERENCE  FUNCTION 


V-PM.V-RM 


X-SCALE-l.OOE-Ql  UNITS  INCH. 
y-SCflLE-2.00E+02  UNITS  INCH. 

CROSS  SPECTRRL  PHASE 


RNGLE,Y-PM,V~RM 


RUN  PH-4 


163 


X-SCALE--=2.00E+01 
Y-SCRLE-2.00E-01 


UNITS  INCH. 
UNITS  INCH. 


TEMPORRL  RUTOCORRELRTION  FN,X-PM,X-RM 


RUN  PH 


4, 


FILE 


( 


OF  C0N6 


164 


oo; 


DC  J 


K-SCRLE-1- OOE 


SCRLE=l-OOE+01 


Gl   UNITS 
UNITS 


INCH, 
INCH. 


POWER   SPECTRUM  -LEUEL    (DB) 
RUN   PH-4,    FILE   7   OF    CGN6 


X-PM,X-RM 


165 


o 

-J 
o 


X-SCRLF----1.00E+00  UNITS    INCH. 
Y-SCRLE-5.00E-02   UNITS    INCH. 

CROSS-CORRELRTIQN   FN, X-PM, X-RM 
RUN   PH-4.    FILE   7   Or    C0N6 


166 


:3H3 


K-SCRLE=i.OOE-Oi  UNITS  INCH 
V-SCflLE=i.OOE+00  UNITS  INCH 

COHERENCE   FUNCTION 


X-PM?X-R 


xH 


X-SCRLE==1.00E-Q1  UNITS  INCH. 
^-SCRL[:--2.00[:+02  UNITS  INCH. 

CROSS  SPECTRRL   PHRSE 


RNGLE,X-PM,X 


PM 
hi  i 


RUN   PH-4 


167 


-K-PM 


I-PY. 


K-oCRLE-l.OOE-Oi  UNITS  IHCH- 
Y-SCRLE-l.OOE+Oi  UNITS  INCH. 

POWER  SPECTRUM  LEUEL  (DB) 
RUN  PH-5,  FILE  8  OF  CPN6 


K-PM,Z-PM 


168 


X-SCRLE-i.OOE+00  UNITS  INCH- 
Y-SCRL£-i .  OOE-Oi  UNITS  INCH- 

CROSS-CORRELATION  FN,X-PM,Z-PM 
RUN  PH-5.  FILE  8  OF  C0N6 


169 


Vx. 


^_ 


-WHfl 


OhJL 


QDd 


X-SCRLE=i.OOE-01   UNITS    INCH- 
V-SCRLE=i.OOE+00   UNITS    INCH. 

COHERENCE   FUNCTION 


X-PM?Z-PM 


^^  \    I —    I  /—v.  I   I  1 — v.  I    I — ^ 


HR? 


I 


X-SCRLE=1.00E-01   UNITS    INCH. 
V-SCRLE-2.00E+02  UNITS    INCH- 

CROSS   SPECTRRL   PHRSE   RNGLE,X-PM 
CH   RLEXRNDER 


"7 

?  L 


PM 


RUN  PH-5 


170 


X-SCRLF-2.00E+01  UNITS  INCH. 
V-SCRLE-2.00E-01  UNITS  INCH. 

TEMPOR.RL  RUTOCQRRELRTION  FN.X-RM 
RUN  PH-5,  FILE  8  OF  CONG 


V 


M 


171 


J?=fiM 


oci 


0C2 


OCo 


UNITS 

UNITS 


K-SCRLE=1-00E-01 
Y-SCRLE=2.00E+Q1 

POWER   SPECTRUM 
RUN   PH-5?    FILE 


IHCM  = 
INCH, 

LEUEL    (DB 
8   OF   C0N6 


K-RM,V-RM 


172 


Q?;Z 


DM 


rt-SCRLF-1.00£>00  UNITS    INCH. 
Y-SCRLE--1.00E-01   UNITS    INCH. 

CRQSS-CCIRRELRT ION   FN ,  X-RM , Y-RM 

RUN   PH-5?    FILE'  8   OF   C0N6 


173 


OOi 


002 


Q23 


004 


JAm 


K-SCRLE=1.00E-Q1   UNITS    INCH 
V-SCRLE=l-OOE+00  UNITS    INCH 

COHERENCE   FUNCTION 

/~\   i  i     /->.  i    i—  ;  i  r~\   i  i  i — ^  i —  i — t 


K-RMpV-RM 


H.qr. 


K-SCRLE=1.00E-01   UNITS    INCH. 
Y-SCRLE--2.00E+02  UNITS    INCH, 


CROSS   SPECTRRL   PHRSE   RNGLE, X-RM, V-RM 


RUN    PH-5 


174 


X-SCALF.=2.00E+0i  UNITS  INCH. 
V-SCflLE-=2.00E-01  UNITS  INCH. 

TEMPORAL  AUTOCORRELATION  FN,V 
RUN  PH-5,  FILE  S  OF  CQN6 


flM,Z-RM 


175 


O 

a 


!  ) 
U 

o 


o:c 


0D2 


004 


X-SCRLE=1.00E-0i   UNITS    INCH- 
Y-SCRLE=1.00E+0I   UNITS    INCH. 

POWER   SPECTRUM   LEUEL    (DB) 
RUN   PH-5,    FILE   8   OF   CQN6 


V-RM.Z-RM 


176 


d 
□ 


x-scfile---\1.ooe+oo  units  inch. 
y-scrle=i.00e-01  units  inch. 

CROSS-CORRELATION   FN,V--RM,Z 

RUN   PH-5,    FILE   8   OF   CONG 


RM 


177 


OCi 


OC2 


iim 


K-5CflLE--=1.00E-01  UNITS  INCH 
Y-SCRLE=1.00E"+OQ  UNITS  INCH 

COHERENCE   FUNCTION 


Y-RN,Z-RM 


HAS 


X-SCRLE=1.00E-Oi  UNITS  INCH. 
y-SCHLE-2.00E+02  UNITS  INCH. 

CROSS   SPECTRAL   PHRSE   RNGLE, Y-RM, Z-RM 

RUN  PH-5 


178 


X-SCRLE=2.0QE+-Q1  UNITS 
V-SCRLE-2.00l-01 


UNITS 


INCH  = 
INCH. 


TEMPORRL   RUTOCORRELRTION   FN,X-RM,Z-RN 
RUN  PH-5,  FILE  8  OF  C0N6 


179 


K-SCRLE-i.OOE-01  UNITS  INCH- 
Y-SCRLE=2.00E+01  UNITS  INCH. 

POWER  SPECTRUM  LEUEL  (DEO 
RUN  PH-5,  PILE  8  OF  C0N6 


x-m,z-m 


180 


X-SCflLE=1.00E+00  UNITS  INCH. 
Y-SCRLE=5. 00E-02  UNITS  INCH. 

CROSS-CORRELRTION  FN, X-RM, Z-RN 
RUN  PH-5,  FILE  8  GF  CQN6 


181 


oj: 


002 


0:3m 


OO-i 


X-SCRLE=1.00E-01  UNJT5  INCH. 
Y-SCRLE=1.00E+00  UNITS  INCH. 

COHERENCE   FUNCTION 


X~RM,Z-RM 


PHRS 


K-SCRLE-l.OOE-01   UNITS    INCH- 
Y-SCflLF=2.00E+02  UNITS   INCH- 

CROSS   SPECTRAL   PHRSE   RNGLE, X-RM, Z-RM 


RUN  PH-5 


182 


:-RM 


K-PH 


K-SCRLE=1.00E-01  UNITS  INCH,. 
Y-SCRLE=i.OOE+Oi  UNITS  INCH- 

POWER  SPECTRUM.  LEUEL  (  DB 
RUN  PH-5p  FILE  8  OF  C0N6 


K-PM?  K-RM 


183 


- 
a 

CI 


a 


X3R 


03C 


032 


X-SCRLE--=i.0OE+00  UNITS  INCH, 
Y-SCflLE=2.OOE-0i  UNITS  INCH. 

CROSS-CORRELATION  FN,X-PM,X-RM 
RUN  PH-5,  FILE  8  OF  C0N6 


184 


X-SCRLE=i.00E-01  UNITS  IHCH 
Y-SCRLE-i.OOE+OQ  UNITS  INCH 

COHERENCE  FUNCTION 


K-PM,X-RH 


PH.qs 


X-SCRLE--=l."OQE-Dl  UNITS' INCH- 
V-SCRLE----2.00E+02  UNITS  INCH, 

CROSS  SPECTRAL  PHRSE  RNGLE,  X-PM-,  X-RM 

RUN  PH-5 


185 


-PK 


X-SCRLE-l.OOE-01  UNITS  INCH. 
V-SCRLE=1.00E+01  UNITS  INCH. 

POUER  SPECTRUM  LEUEL  ( DB 
RUN  PH-5,  FILE  8  OF  CQN6 


Z-PM.Z-RM 


186 


O 


a 
o 


—i 

o 
I 


do-; 


CCGR 


K-SCHLE--1.00E+00  UNITS  INCH. 
V-SCRLE=i.00E-01  UNITS  INCH. 

CROSS-CORRELRTION  fn,?-pm,z-rm 

RUN  PH-5,  FILE  8  OF  C0N6 


187 


:wr 


021 


322 


QQ1 


02-1 


X-SCRLE=J..0QE-Q1  UNITS  INCH 
V-SCRLE=1.00'E+00  UNITS  INCH 

COHERENCE   FUNCTION 


Z-PM,Z-RM 


vPHSS 


X-SCflLE=l.QQE-Qi  UNITS   INCH 
V-SCRLE-2.00E+02   UNITS    INCH 


CROSS   SPEGTRRL   PHRSE   RNGLE, Z-PM, Z-flM 


f~\    I   I      S— v.  I 


RUN  PH-5 


188 


X-SCRLE=2.00E+0i   UNITS    INCH- 
Y-SCRLE=2.00E-0i   UNITS    INCH. 

TEMPORRL   RUTQCORRELRTION   FN, V-PM, Z-PM 

RUN   PH-6,    FILE   9   OF   CQN6 

1~   S.iiV"* 


189 


X-SCRLE=i.OOE-Oi  UNITS  INCH. 
V-SCRLE-l.OOc+Ol  UNITS  INCH. 

POWER  SPECTRUM  LEUEL  (DB) 
RUN  PH-6,  FILE  9  OF  C0N6 


V-PM.Z-PM 


190 


o:c 


OCi 


0:2 


rim 


0D4 


K-SCRLE=1.00E+00 
y-SCRLE=2.00E-01 

CROSS 

RUN   PH-6,    FILE 


UNITS    INCH. 
UNITS    INCH- 

FN,Y-PM,Z 
]F   C0N6 


CQRRELRTION 


PM 


191 


X-SCRLE=1. OOE-01  UNITS  INCH 
Y-SCRLE-LOOE+OQ  UNITS  INCH 

COHERENCE  FUNCTION 


Y~PM,Z-PM 


K-5CRLE=1.O0E-O1  UNITS  INCH. 
Y-SCRLE=2-QQE+0'2  UNITS  INCH. 

CROSS  SPECTRRL  PHRSE  RNGLE, V-PM, Z-PM 

RUN   PH-6 


192 


\S     V  -PM 


302 


001 


006 


009 


X-SCflLE=2.00E+0I  UNITS  INCH. 
y-SCfiLE=2.00E-01  UNITS  INCH. 

TEMPORAL  nUTOCQR'RELflTION  FN.X-PM.Z-P 

RUN  PH-6.  FILE  9  OF  C0N6 


193 


K-Pl-1 


S.~PX 


K-SCnLE=1.00E-0J   UNITS    INCH. 
Y-SCRLE-5-  OOE+00  UNITS.   INCH. 

POHER   SPECTRUM  LEUEL    (DB) 
RUN  PH-6,    FILE   9   OF    C0N6 


X-PM,Z-PM 


194 


,fcC3R 


o:c 


032 


X-SCRLE-LOOE+OO  UNITS    INCH. 
Y-SCRLE=1.00E-01   UNITS    INCH. 

CROSS-CORRELRT I  ON   FN , K-PM, Z-PM 
RUN   PH-6,    FILE   9   OF   C0N6 


195 


_LDHT 


X-SCRLE=1.00E-01  UN1T5  INCH 
V-SCRLE=l-OOE+00  UNITS  INCH 

COHERENCE  FUNCTION 


K-PM.Z-PM 


X-SCflLE=l.QOE-01  UNITS  INCH. 
Y-SCRLE=--2-00E+02  UNITS  INCH> 

CROSS  SPECTRRL  PHRSE  RNGLE, X~PM, 2-PM 

RUN  PH-6 


196 


X-SCRLE=--2.00E+01  UNITS  INCH- 
y-SCRLE--=2.00E-01  UNITS  INCH. 

TEMPORAL  RUTQCORRELRTIQN  FN, X-PM, X-R 
RUN  PH-6,  FILE  9  OF  C0N6 


197 


k-ph 


x-scfile=i.  ooe-oi  units  inch. 
y-scfile--5.  ooe+00  units  inch. 

POWER   SPECTRUM  LEUEL    (DB) 

RUN   PH-6,    FILE   9   OF   C0N6 


X-PM,X-RM 


198 


oo: 


Oj-1 


:zn 


K-SCRLE==1.00E+00  UNITS  INCH. 
Y-SCRLE=5.00E-02  UNITS  INCH- 

CROSS-CORRELATION   FN,X-PM,X-RM 
RUN   PH-6,    FILE   9   OF   C0N6 


199 


vLGH'i 


K-SCRLE=1.0QE-Q1   UNITS    INCH. 
V-SCRLE=l-OOE+00   UNITS    INCH. 

COHERENCE   FUNCTION      K-PM.K-RM 


X-SCRLE=1.QQE-01   UNITS    INCH- 
^-SCnLE----2.00E:+02   UNITS    INCH. 

CROSS   SPECTRAL   PHRSE   RNGLE, X-PM,X-RM 

RUN    PH-6 


200 


y-PM 


K-SCRLE=2.00E+01 
Y-SCRLE-2- OOE-01 


UNITS  INCH. 
UNITS  INCH. 


TEMPORAL   AUTOCORRELATION   FNpV-PM,Y-AM 

RUN   PH-6?    FILE   9   OF   C0N6 


201 


V-n:-: 


t>  __  '  _ 


V.pw 


K-SCflLE-a.OOE-OJ;  UNITS  INCH. 
V-SCRLE-l-OOE'+Ol  UNITS  INCH. 

POWER   SPECTRUM   LEUEL    (DB) 
RUN   PH-6,    FILE   9   OF   C0N6 


Y-PM,Y-RM 


202 


K-SCRLE--1.  OOF  +  00  UNITS  INCH- 
V-SCRLE==i.OOE-Gi  UNITS  INCH. 

CROSS-CORRELATION  FN,Y-PM,V-RM 
RUN  PH-6,  FILE  9  OF  C0N6 


203 


^-/X. 


QZA 


X-SCRLE==1.00E-Q1  UNITS  INCH. 
Y-SCRLE=l.Q0E+00  UNITS  INCH- 

COHERENCE  FUNCTION  Y~PM,Y~RM 


HI  3 


K-5CRLE=1-0QE-01  UNITS  INCH- 
Y-SCRLE=2.00E+02  UNITS  INCH- 

CROSS   SPECTRAL    PHRSE   RNGLE,  V-PM,  SJ-RM 

RUN    PH-6 


204 


V-PM 


•;-.rjM 


034 


006 


038 


X-SCALE----2.00E+01  UNITS  INCH- 
Y-SCflLE=2.00E-01  UNITS  INCH. 

TEMPORRL   AUTOCORRELATION   FN,X-PM,Y-PM 

RUN   PH-6,    FILL7   9   OF   C0N6 


205 


V-PM 


K-5CRLE=i. OOE-01  UNITS  INCH. 
V-SCRLE--5.00F  +  00  UNITS  INCH. 

POWER  SPECTRUM  LEUEL  (DB) 
RUN  PH-6,  FILE  9  OF  C0H6 


K-PM-V-PM 


206 


L'VW 


002 


X-SCflLE-LOOE+OO  UNITS  INCH. 
Y-SCHLE--=1.00E-01  UNITS  INCH- 

CROSS-CORRELATION  FN,X-PM,Y-PM 
RUN  PH-6,  FILE  9  OF  C0N6 


207 


UC2 


034 


c<2>JSj^ 


K-SCRLE=1.D0E-01   UNITS    INCH 
Y-SCRLE=1.00E+Q0  UNITS    INCH 

COHERENCE   FUNCTION 


X-PM,Y-PM 


o 
o 


vHRS 


X-SCRLE=1.00E~Qi  UNITS   INCH- 
Y-SCRLE=--2-00E+02  UNITS   INCH. 


CROSS   SPECTRAL   PHRSE   ANGLE, X-PM, V-PM 

RUN    PH-6 


208 


'KM" I  I    J  ' 


r 


UCS 


K-i C-LI  -2-00E+01  UNITS  INCH- 
Y-SCRLF.--=2.00E-01  UNITS  INCH. 

<TEMPORRL   nUTOCORRELRTIGN   FN,X-PM   V-Ph 
RUN   PH-7,    FILE   10   OF   CON   6 


209 


K~SCRLE----1.00E-0i  UNITS  INCH. 
V-SCRLF.----J..00E+01  UNITS  INCH. 

POWER  SPECTRUM  LEUEL  (DB)  K-PM  V-Ph 
RUN  PH-7,  FILE  10  OF  CON  6 


210 


a 
a 
a 


00<l 


X-SCRLE^LOOE+OO  UNITS    INCH- 
Y-SCRLE=--1.00E-0i   UNITS    INCH. 

CROSS-CDRRELRTIGN   FN, X-PM, V-D 
RUN   PH-7,    FILE    10   OF   CON   6 


?  i 


211 


V 


r 


002 


004 


K-SCRLE-LOQE-Ql   UNITS   INCH- 
Y-SCRLE=1.0QE+00   UNITS    INCH. 

COHERENCE   FN,    X-PM?Y-PM 


•;H5 


K~SCRLtr-l  =  GOL:.-Gi  UNITS  INCH 
V-SCRLE=2.Q0E+02  UNITS  IN 


L»ri« 


CROSS   SPECTRAL   PHRSE   RNGLE,X-PM   ^'-Pi- 


rn   PH-7 


212 


>-G(  Hi  C^2.GGF+01   UNITS    INCH- 
Y-SCRLE-2-00E-01  UN  J  fS    INCH. 


TEMPORRL  RUTOCORRELRTION  FN,V-PM,V-R! 


■- 


RUN 


Ml 


f   9 


LE  10  Or  CON  6 


213 


y-PM 


?,r>. 


ICH- 


*• 


-i.OOF-01  UNITS  I 
SCRLE=-i.  OOE+Oi  UNITS  INCH. 

PQNER  spectrum  leuel 

RUN  PH-7,  FILE  10  OF 


(DB)  Y-PM,V-RM 

CON  6 


214 


a 

a 


X-SCRLE=  L.OOE+00  UNITS    IN! 
V-SCRLE=1.00E-01   UNITS    INCH. 

CROSS-CORRELRTION   FN,Y-PM,Y 

RUN  PH-7.    FILE   10   Or    CON   6 


Rh 


215 


I  1 

:  i 


\       tJ 


./ 


vA 


oj: 


D02 


00'5 


X-SCRLE=1.00E-01   UNITS    INCH, 
V-SCRLE---1.00E +00  UNITS    INCH. 

COHERENCE   FUNCTION 


Y-PM/Y-RM 


□ 

a 
a 


PHflS 


X-SCRLE-1.09E-01   UNITS    INCH- 
Y-SCRLE=2.Q0E+02  UNITS   INCH, 

CROSS   SPECTRAL   PHRS 


RNGLE.Y-P 


)M     \, 


-v 


RUN    PH-7 


216 


J<-PM 


K-SCRLE==2.00E+01  UNITS  INCH, 
Y-SCRLE=2-0QE-01  UNITS  INCH. 

TEMPGRRL  RUTQCQRRELRTION  FN.K-RH 
RUN  PH-7?  FILE  10  OF  CON  6 


Y-fi 


M 


217 


«•  -RM 


000 


001 


002 


003 


::< 


X-SCRLE=i.00E-01   UNITS    INCH. 
Y-SCRLE=2.00E+0i   UNITS    INCH. 

PQNER   SPECTRUM   LEUEL 


RUN   PH 


"1 

f  9 


FILE' 10   OF 


(DB)    X-RM,Y-R 


M 


CON   6 


218 


( > 
rn 
D 


CQR 


0:2 


001 


X-SCRLE-l.OOE+0'0  UNITS  INCH- 
Y-SCflLE=--i.OOE-Dl  UNITS  INCH. 

CRGSS-CORRELRTION  P1 

RUN  PH-7,  FILE  10  OF 


X-RM, 

CON  6 


y 


RM 


219 


x_ 


GCi 


yv 


0Q2 


' 


X-SCRLE=1.00E-0J  UNITS  ]  H^:  i, 
Y-SCRLE= l.QOE+00  UN3  rs  INCH. 

COHERENCE  FUNC1 ION 


X-RMpV-  RM 


X-SCALE-1.00E-01  UNITS  INCH. 
Y-SCRLE-2.00E+Q2  UNITS  INCH. 

CROSS  SPECTRRL  PHRSE  RNGLE,  X-RM',  V-F 

RUN   PH-7 


220 


/>-^1 


K-SCRLE=I.00E-01  UNITS  INCH. 
V-SCRLE^l.OOE+01  UNITS  INCH. 

POhlR  SPECTRUM  1EUEL  I 03 
RUN  PH-7.  FILE  10  Of:  CON 


X-PM 


n 


o 


221 


X-SCRLE>i.OUE+00  UNITS  INCH. 
Y-SCRLE----5.00E-02  UNITS  INCH. 

CRQSS-CQRRELRTICjN  FN,X-PM,X 
RUN  PH-7,  FILE  10  OF  CON  6 


RN 


222 


ac ; 


002 


QL  3 


x      ' 


K-SCRLE-1.  00E-01  UNITS  INCH- 

v-scrle-i.ooe+qo  units  inch. 

COHERENCE   FN.    X-PM,X-RM 


K-SCRLE-i-OQE-Oi  UNITS  INCH. 
V-SCRLE=2-00E+02  UNITS  INCH. 

CROSS  SPECTRRL  PHRSE 


flNGLE,X-PM,X-RM 


RUN   PH-7 


223 


oco 


002 


UC"'i 


CC3 


X-SCflLE=2.00E+0i  UNITS  INCH. 
Y-SCRLE-2.00E-01  UNITS  INCH. 


TEMPORAL   RUTOCORRELRTION   FN.V- 
RUN   PH-3      FILE    11   Op   CON   6 


224 


K-SCRLE----1.00E-Q4.  UNITS    INCH. 
Y-SCRLE---2-00E+01   UNITS    INCH- 

POWER   SPECTRUM   LEUEL 
RUN   PH-8,    FILE   11   VL 


(DB) 

CON   6 


'v1 


p 


"7 


PM 


225 


nnr 
UO*. 


aci 


002 


003 


004 


X-SCRLE=1.QQE+Q0  UNITS    INCH. 
Y-SCRLE=2-00E-01   UNITS    INCH- 

CROSS-CQRRELRTIGN   FN,V-PM, 
RUN   PH-E    FILE   11   3^   CON   5 


i 


226 


0J2 


D3«j 


X-SCRLE=1.QQE-01  UNITS  INCH 
Y-SCRLE--=l-QOF+00  UNITS  INCH 

COHERENCE  FUNCTION 


Y-PM  RND  Z-PM 


K-SCRLE=l-OOE-Ql  UNITS  INCH. 
V-SCflLE--=2.00E+02  UNITS  INCH. 

CROSS  SPECTRRL  PHRSE  ANGLES 


— 


RUN  PH-8 


227 


K-SCRLF-l.QOE-01  UNITS  INCH. 
y-S'.^Lt---1  . OOE+01  UNITS  INCH. 

POWER   SPECTRUM'  LEUEL    (DB) 
RUN  PH-8,  FILE  11  OF  CON  6 


PM  } 


228 


K-SCRLE----1-  OQr+OO  UNITS  INCH. 
y-SCRLE--=2.00E-01  UNITS  INCH- 

CROSS-CORRELnTlON  FN, Z-PM, Z-RM 
RUN  PH-11  FILE  11  OF  CON  6 


229 


a  V  / 

K  i 


OCJ 


r  02 


V 


X/\z 


■ 


K-SC  !LE=1. 00E-Q1  UNITS  INCH 
Y-SCRLE==1.0QE+00  UNITS  INCH 

COHERENCE   FUNCTION 


Z-Pil   RMD 


L 


K-SCflLE--=1.00E-01  UNITS,  INCH. 
y-SCflLE-2.00E.-f02  UNITS  INCH. 

CROSS  SPECTRRL  PHRSE 


RNGLE,Z-PM,Z-RM 


RUN   PH-8 


230 


u 


■ 


~V-PM 


■  :  :RLE--=i.  OOE-Oi  UNITS  INCH. 
V-SCRLE--2.Q0E+01  UNI  r5  INCH- 

POWER  SPECTRUM  LEUEL  ([ 
RUN  PH-8,  EIL::  11  OF  CON 


^ 


'-PM.S 


N  . 


o 


231 


,C~P 


DZZ 


d:i 


K-5CRLE-1.00E+00  UNITS  INCH. 
Y-SC.RLE:=1.00E-0i  UNITS  INCH- 

CROSS-CORRELRTION   FN,Y-PM,V-RM 
RUN   PH-B    FILE   11   OF   CON   6 


232 


iL 


£^ 


GOZ 


OC4 


-  Vi 


X-SCRLE-l.OOE-Ol  UHJTS  INCH. 
V-SCRLE--1.00E+00  UNI  IS  INCH, 

COHERENCE  FUNCTION 


v  -  i 


• RND  V-RM 


K-SCRLE----1-OOE-01  UNITS  INCH. 
Y-SCRLE-2.00E+Q2  UNITS  INCH- 

CROSS  SPECTRRL  PHASE  RNGLE.V-PM  V-RM 

RUN  PH-8 
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